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profile  is  from  a  fit  to  SIMS  data  6 
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3.2.  (a)  Dopant  concentrations  and  abrupt-depletion  electron 

profiles  vs  depth.  The  corresponding  energy-band 
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EXECUTIVE  SUMMARY 


I  . 

The  Advanced  Growth  and  Process  Modeling  for  GaAs  Program  was 
conducted  by  the  GaAs  process  modeling  and  technology  group  of  Professor 
T.W.  Sigmon  in  the  Solid  State  Laboratory  of  the  Electrical  Engineering 
Department  at  Stanford  University.  The  period  of  performance  began  in 
December  1985  and  ended  in  October  1988.  This  document  is  the  final  report  of 
the  program  and  presents  and  discusses  the  major  achievements  and 
significant  findings. 

The  objective  of  the  program  was  to  develop  physical  models  for 
substrate  and  process  defects  and  variations  which  affect  the  electrical 
parameters  of  MESFET  structures  fabricated  on  this  material  by  direct  ion 

implantation  processes.  Verification  of  these  models  was  to  be  carried  out 
where  possible  by  fabrication  of  specially  designed  device  structures  capable 
of  providing  statistically  significant  data.  In  order  to  remove,  as  much  as 
possible,  unwanted  process  variations  all  device  processing  was  carried  out 
under  a  subcontract  with  the  Rockwell  MRDC  GaAs  pilot  line  facility. 

The  technical  approach  of  the  program  was  to  concentrate  on  direct  ion 
implanted  MESFET  process  technology  fabricated  on  EL2  compensated  semi- 
insulating  substrates.  This  decision  was  supported  by  the  improved 

manufacturbility  provided  by  this  approach  over  epi-growth  processes. 
Following  our  acessment  of  this  technological  approach  to  fabricate  MESFET 
circuits  and  to  complement  other  programs  being  administered  by 
AFWAL/MLPO  the  following  specific  areas  were  singled  out  for  investigation. 

1. )  Development  of  models  to  determine  the  effect  of  substrate 

impurities  and  deep  levels  (EL2,  Cr)  on  MESFET  electrical 

parameters. 

2. )  Development  of  analytical  models  to  calculate  ion  implantation 

profiles  in  GaAs  including  both  channeling  and  encapsulation 

effects. 

3. )  Development  of  a  process  model  (GATES)  to  allow  calculation  of 

the  effects  and  study  the  sensitivity  of  the  above  effects  on 
MESFET  electrical  parameters. 

4. )  Development  of  a  physical  model  for  how  dislocations  affect 

MESFET  electrical  parameters,  extension  of  this  model  to  wafer 
scale  by  the  use  of  Monte  Carlo  techniques. 


5.)  Design  of  a  D_ense  R_ow  P.attern  (DRP)  of  MESFETs  and  the 

appropriate  automatic  test  algorithms  for  verification  of  the 
predictions  of  the  above  models. 

The  primary  successes  of  the  program  include:  i.)  development  of  a 
physical  model  to  explain  the  effect  of  dislocations  on  the  fluctuations  in 
threshold  voltage  observed  experimentally,  ii.)  development  of  analytical 
Pearson  IV  models  for  the  profiles  of  the  technologically  significant  ions  Si, 

Se,  and  Be  in  <100>  oriented  GaAs  for  specific  tilt  and  rotation  angles  of  the 
substrates  with  respect  to  the  implant  beam,  iii.)  development  of  a  one 
dimensional  device  and  process  simulator  which  allows  exploration  of  the 

effect  of  certain  process  variables  on  the  significant  MESFET  electrical 

parameters. 
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1 1 .  MECHANISM  OF  EL2  EFFECTS  ON  GaAs  FIELD-EFFECT  TRANSISTOR 

THRESHOLD  VOLTAGES 

2.1  Background 

In  recent  years  much  effort  has  been  devoted  to  characterizing  the 
concentration  of  the  deep-donor  EL2,  present  in  semi-insulating  liquid- 
encapsulated-Czochralski  (LEC)-grown  GaAs  (1-3).  The  average  EL2 
concentration  varies  across  wafers,  possibly  originating  from  changes  in  the 
thermal  stress  during  boule  growth  (4).  Also,  local  fluctuations  in  the  EL2 
concentration  over  50-  to  SOO-pm  distances  have  been  observed  around 
individual  dislocations  and  around  clusters  or  lines  of  dislocations  (1-3). 

These  fluctuations  in  EL2  concentrations  are  a  possible  candidate  for 

explaining  threshold  voltage  fluctuations  in  ion-implanted  metal- 

semiconductor  field-effect  transistors  (MESFETs).  Dobrilla  et  al.  (2)  measured 
threshold  voltages  in  ion-implanted  MESFETs  where  the  EL2  concentration 
varied  from  7  to  10x10^^  cm'3.  A  linear  correlation  was  found;  Vth  decreased 
by  63  mV  per  lO^^  EL2/cm^,  though  the  magnitude  of  the  variations  over  the 
range  of  3xl0l5cm’3  were  as  small  as  the  random  fluctuations  in  threshold 
voltages. 

In  this  work  we  discuss  a  mechanism  for  EL2  effects  on  threshold 
voltages  which  focuses  on  its  role  in  fixing  the  substrate  barrier  height.  We 

are  not  so  much  concerned  with  the  presence  of  EL2,  but  with  its  absence,  due 
to  outdiffusion. 

2.2  Model  and  Calculations 

Makram-Ebeid  et  al.  (6)  found  that  for  capless  annealing  with  ASH3 

overpressures  at  temperatures  between  750  and  900°C,  the  EL2  concentrations 
varies  as 

[EL2]  =  [EL2I00  erf  ^x/2v/Dt  j,  ^2.1) 

where  2^^  is  between  0.5  to  5  pm  for  20  min  anneals.  F  or  other  annealing 
caps,  similiar  diffusion  lengths,  but  more  complicated  profile  types  were 

observed.  Since  most  MESFETs  are  fabricated  with  implant  depths  x<l  pm,  and 
arc  annealed  for  10-30  min  at  these  temperatures,  EL2  outdiffusion  can  be 

present  to  some  degree  in  all  furnace-annealed,  ion-implanted  MESFETs. 
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To  explain  our  model  we  consider  the  energy  band  diagrams  associated  with 
two  possible  EL2  outdiffusion  profiles,  EL2  and  EL2',  in  Fig.  2.1.  For  EL2', 
annealing  for  a  longer  time  or  higher  temperature  is  assumed.  At  the 


Fig.  2.1.  (a)  Concentration  profiles  of  implanted  donor  (Nd  ) . 

background  acceptor  (Na).  and  two  different  EL2  profiles. 
The  shaded  areas  show  the  substrate  and  Schottky  barrier 
depletion  regions.  (b).  Energy  band  diagrams  for  the 
profiles  EL2  and  EL2’  in  1(a). 

intersection  xq,  between  the  implanted  shallow  donor  and  the  background 
shallow  acceptor  [usually  carbon,  present  in  concentrations  of  1  to  6x10^  ^ 
cm  *3  in  LEC  GaAs  (7)],  a  depletion  region  forms  around  a  barrier  of  height 
Vs-, 5.  The  magnitude  of  Vsub  depends  on  the  differences  between  A2=Ec-Ef  on 
the  substrate  side  of  the  region  and  A 1  (~0.05  V)  on  the  channel  side.  For  EL2 
concentrations  higher  than  the  acceptor  concentration  on  the  substrate  side. 
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the  Fermi  level  is  pinned  to  the  mid-gap  EL2  energy  level,  so  that  A 2-0.75  V  and 
Vsub~0.7  V.  When  EL2’  is  less  that  the  the  acceptor  concentration  on  the 
substrate  side,  the  Fermi  level  will  be  closer  to  the  valence  band  so  that 
A2~^sub  =  l-2  V~Egap.  The  larger  substrate  barrier  height  in  the  latter  case 
pushes  the  depletion  region  in  the  n-Iayer  toward  the  surface,  reducing  the 
width  of  the  channel,  and  resulting  in  less  negative  threshold  voltages  being 
needed  to  fully  deplete  the  channel. 

Therefore  there  are  two  possible  ways  EL2  can  affect  threshold  voltages. 
If  EL2  is  initially  dominant  on  the  substrate  side,  annealing  for  longer  times 
can  shift  the  substrate  barrier  height,  and  therefore  Vth-  If  the  anneal  time 
at  the  average  bulk  EL2  concentration  is  such  to  give  a  substrate  barrier 
height  between  0.7  and  1.2  V,  local  fluctuations  in  the  bulk  EL2  concentration 
can  give  fluctuating  threshold  voltages. 

To  test  these  ideas,  two  model  calculations  are  made.  In  the  abrupt- 

depletion  model,  (5)  we  integrate  over  the  substrate  depletion  region  and  do  a 
three-component  Fermi-level  anaysis  on  the  substrate  side  to  determine  self 
consistently  how  large  Vgub  ihe  depletion  widths  must  be  to  satisfy 

Vsub  =  ^2'^l-  Ifi  model,  the  contribution  of  the  ionized  EL2  in  the  depletion 
region  is  neglected;  only  its  effect  in  determining  the  Fermi-level  position  at 

the  edge  of  the  depletion  region  is  included. 

In  a  second  calculation,  Poissons  Equation  (8)  is  solved  for  the  electron 
concentration  n(x),  including  the  shallow  donor  and  acceptor  and  partially 
ionized  EL2.  No  bias  is  applied,  and  as  a  boundary  condition,  dn(x)/dx  is 
assumed  to  vanish  at  the  surface.  Figure  2.2  shows  the  calculated  Nd(x),  Na(x), 
and  electron  concentrations  n(x)  for  various  EL2  outdiffusion  profiles.  The 
shape  of  the  implanted  donor  profiles  is  not  Gaussian,  but  has  a  significant 
channeling  tail  (Sec  Section  IV)  measured  using  secondary-ion  mass 
spectroscopy  (SIMS)  (9).  For  Si  ions  implanted  with  energies  less  than  200 

keV,  these  tails  are  nearly  independent  of  the  wafer  orientation  for  tilt  angles 

greater  than  4°.  However,  the  electron  concentration  profiles  are  nearly 
always  Gaussian.  In  the  abrupt-depletion  approximation,  n(x)  vanishes  for 
x>xi,  where  xt  is  the  channel  edge  of  the  substrate  depletion  region.  In  the 
Poisson  solution,  n(x)  does  not  vanish,  but  shows  a  Debye  tail.  In  the 
literature,  LSS  donor  profiles  are  often  compared  with  electron  profiles, 
measured  using  the  capacitance-voltage  (CV)  technique  (10).  Figure  2.2 
demonstrates  that  the  apparent  good  agreement  between  Gaussian  LSS  implant 
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Fig.  2.2.  Donor  Ndi  acceptor  Na.  and  calculated  electron  concentration 
profiles  for  several  different  outdiffused  EL2  profiles,  assuming 
implanted  60  keV  Si  ions  with  a  dose  of  2x10*  ^  ions/cm^.  (The 
anneal  time  in  minutes  and  the  bulk  EL2  concentration  in 
10*®  cm *3  arc  given).  The  positions  of  the  substrate  junction  Xq 
and  the  shallow  edge  xt  of  the  substrate  junction  in  the  abrupt 
depletion  approximation  are  shown  by  arrows.  The  implanted 
profile  is  from  a  fit  to  SIMS  data. 


profiles  and  CV  profiles  occurs  because  ion  implantation  into  EL2  (or  Cr) 
compensated  GaAs  is  like  implantation  into  p-layers;  the  donor  channeling 
tails  are  consumed  by  the  background  acceptor  concentrations. 

For  the  EL2  profiles  used  to  calculate  the  electron  profiles  shown  in  Fig. 
2.2,  we  assumed  an  error-function-type  outdiffusion  with  a  diffusion  constant 
D  =  3xl0'**  cm^/sec,  measured  for  an  800°C  capless  anneal  with  AsH3 
overpressure  (6).  As  the  annealing  time  increases,  Vgub  increases  from  0.7  to 
1.2V,  and  the  electron  concentration  profiles  are  narrower. 

Figure  2.3  shows  the  EL2  effects  on  threshold  voltages.  Using  the  Poisson 
solution  of  n(x),  Vth  is  calculated  using 


_5_ 


xdx  n(x), 

0 


(2.2) 
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where  Vbi  is  the  Schottky  barrier  height.  With  the  abrupt-depletion 
approximation,  one  must  first  calculate  xt  by  integrating  over  the  substrate 
barrier.  Then  Vth  is  obtained  from  (5) 


Vth  =  Vbi 


pXt 


xdx  (nND(x)-NA.(x))  -  A^, 


(2.3) 


where  is  the  donor  activation  (assumed  unity  in  Fig.  2.3). 


^  EL2(10^®cm*^) 

(b) 


Fig.  2.3.  (a)  Calculated  threshold  voltages  for  60  keV  Si  ions  vs  anneal 
time  for  a  constant  bulk  EL2  and  acceptor  concentration  of 
2x]0'^and  4x10* ^  cm'3,  respectively.  (b)  Vth  vs  bulk  EL2 
concentration  for  a  20  min  anneal  at  800°C. 
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In  Fig.  2.3a,  as  the  annealing  time  tA  increases,  Vth  increases  rapidly  in  the 
abrupt-depletion  model  near  tA  ~  20  min.  This  is  the  point  where  the  EL2 
concentration  falls  below  Na  on  the  deep  side  of  the  substrate  depletion 

region.  For  tA<10  min  and  tA>30  min.  Vgub  is  either  0.7  or  1.2  V,  and  V^h  does 
not  vary  much  with  annealing  time.  With  the  Poisson  solution,  the  change  is 
less  rapid,  and  occurs  over  the  whole  range  of  annealing  times,  due  in  part  to 
the  inclusion  of  the  partially  ionized  EL2  in  determining  the  depth  of  the 
substrate  depletion  region.  If  tA  is  constant,  and  the  EL2  bulk  concentration 
increases,  Vth  decreases  with  both  solvers,  as  shown  in  Fig.  2.3b.  When  the 
bulk  EL2  concentration  falls  below  10^^  cm'3,  Vth  increases  by  about  75  mV. 

2.3  Conclusion 

We  conclude  that  for  special  cases,  local  fluctuations  in  EL2 

concentrations  from  3  to  20x10^^  cm '3  could  possibly  affect  MESFET  threshold 
voltages.  However,  the  effect  is  never  sufficiently  large  to  explain  the 

measurements  of  Dobrilla  et  al.  (2).  The  largest  effect  we  calculated  for  the 
MESFETs  measured  (Vth=-200  mV)  is  a  shift  of  13  mV/10" EL2/cm3,  compared 
with  63  mV  measured.  The  calculations  assumed  a  15  min  anneal  at  800°C, 
error  function-type  outdiffusion,  D=3xl0"^*  cm^/sec,  implanted  60-kcV  Si  ions 
with  a  dose  of  1.5x10*2  ions/cm^,  =90%,  and  Na=4x10*5  cm'3.  The  magnitude  of 
the  effect  depends  sensitively  on  the  anneal  time;  the  shift  is  less  than  5 
mV/1 0*5  EL2-cm"3  fortA<10  min  and  tA>25  min.  If  this  were  the  major 
mechanism  for  threshold-voltage  fluctuations  in  MESFETs,  it  would  be  fairly 
easy  to  eliminate  by  preannealing  the  wafers  before  implantation  for  long 

times  to  drive  out  the  EL2,  using  a  buried  p-layer  (11)  to  fix  Vjub  to  ~1.2  V,  or 

by  annealing  the  whole  boule  for  a  long  time  to  obtain  a  uniform  EL2 
concentration  (12). 

By  comparing  the  calculated  magnitude  of  EL2  effects  with  reported 
threshold-voltage  standard  nonuniformities  (13),  we  conclude  that  this 
mechanism  is  probably  not  the  most  important  mechanism  for  giving  MESFET 
threshold-voltage  fluctuations;  other  mechanisms  such  as  the  dislocation- 
increased  carrier-activation  (DICA)  model  (13)  are  less  fragile  and  give  larger 
fluctuations. 
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III.  SUBSTRATE-IMPURITIES  EFFECTS  ON  GaAs  MESFETs 


3.1  Background 

Materials  nonuniformities  in  Liquid-Encapsulated-Czochralski  (LEC) 
grown  GaAs  substrates  have  caused  difficulties  in  the  GaAs  semiconductor 
industry.  There  are  two  important  classes  of  materials  effects;  the  effect  of 
dislocations  [13-15]  and  of  background  impurities  [16-19].  We  discussed 
calculations  of  the  effects  of  dislocations  on  the  properties  of  ion-implanted 
GaAs  MESFETs,  in  our  Imterim  Report. 

While  dislocations  can  have  a  local  effect  on  MESFET  threshold  voltages 
and  other  device  characteristics,  impurities  will  affect  device  characteristics 
approximately  uniformly  across  an  entire  wafer  [16]  when  they  are 
uniformily  distributed.  As  typical  GaAs  MESFETs  are  fabricated  with  active 
channel  donor-impurity  peak  concentrations  of  about  lO'^  cm'^,  background 
impurity  levels  above  10^^  cm*3  may  affect  device  characteristics.  In  LEC- 
grown  GaAs,  shallow  acceptor  levels  of  carbon  are  present  which  vary  from  1 
to  6  X  10^^  cm'3  [18].  Background  donor  levels  of  S  or  Si  are  also  present,  but  at 
concentrations  well  below  10^5  cm'^,  and  usually  below  the  C  concentration 
[17,19].  Bulk  LEC  grown  GaAs  is  semi-insulating  due  to  tlte  presence  of  the 
deep-donor  EL2  at  levels  of  about  10^^  cm"3  [19].  Excess  EL2  over  the  shallow 
acceptor  pins  the  Fermi  level  at  the  middle  of  the  energy  gap,  resulting  in 
small  electron  concentrations  and  high  resistivity  in  the  intrinsic  material. 
In  Horizontal-Bridgman  (HB)  grown  GaAs,  an  excess  of  background  shallow 
donors  over  shallow  acceptors  is  usually  present,  and  the  deep  acceptor  Cr  is 
added  to  pin  the  Fermi  level  to  midgap.  Unfortunately,  both  Cr  [20]  and  EL2  [6] 
outdiffuse  in  the  near  surface  at  temperatures  characteristic  of  the  post¬ 
implant  anneals.  For  this  reason,  Cr-doped  wafers  have  been  abandonned  by 
many  of  the  MESFET  manufactures  in  favor  of  the  LEC-grown,  EL2-pinned  GaAs 
substrates.  The  outdiffusion  of  the  EL2  can  have  an  effect  on  GaAs  MESFET 
threshold  voltages  [21],  but  it  is  relatively  minor  when  compared  with  that  of 
Cr  outdiffusion. 

The  incorporation  of  impurities  from  the  GaAs  melt  generally  varies 
from  wafer  to  wafer,  and  from  boule-to-boule  so  that  impurity  effects  can 
result  in  boule-scale  nonuniformities  in  device  characteristics.  To  reduce 
these  effects  it  is  common  to  qualify  boules.  Test  MESFET  or  Schottky  diode 
structures  are  fabricated  on  wafers  taken  from  near  the  seed  and  tail, 
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threshold  voltages  are  measured,  and  the  ion-implant  doses  needed  to  obtain 

uniform  device  characteristics  are  calculated  [22],  Since  the  segregation  of 
impurities  varies  regularly  along  the  growth  direction  [16],  the  doses  needed 
to  obtain  the  desired  device  characteristics  on  the  remaining  untested  wafers 
of  the  boule  are  simply  interpolated. 

It  is  clearly  desirable  to  reduce  impurity  levels  as  much  as  possible,  to 
below  10^^  cm‘^,  to  minimize  impurity  effects  on  MESFET  devices.  However, 

eliminating  background  acceptor  impurities  in  LEC  GaAs  is  undesirable,  since 
these  background  acceptors  minimize  the  channeling  tails  of  the  implanted 

donor  ion  distributions.  If  the  wafer-to-wafer  impurity  levels  cannot  be 
tightly  controlled,  the  next  best  thing  to  do  is  to  eliminate  background 

impurities  altogether.  However,  if  the  wafer  manufacturers  succeed  in  this,  it 
will  be  necessary  for  the  device  manufacturers  to  add  a  buried  p  layer  to 
obtain  sharp  electron  profiles,  low  substrate  leakage,  and  device  isolation. 

In  Section  3.2  we  discuss  the  results  of  model  calculations  of  impurity 
effects  on  GaAs  MESFET  threshold  voltages  and  electron  profiles.  The 

dependence  of  threshold  voltages  on  acceptor  concentrations  is  found  to  be 
significant.  To  reduce  this  dependence,  buried  p  layers  (Sect.  3.2),  fabricated 
by  implanting  Be  ions  can  be  used.  Calculations  are  presented  which  vary  the 

Be  ion  implant  energy  and  dose,  and  calculate  the  effect  of  varying 

background  C  concentrations.  In  Section.3.3  we  discuss  the  conclusions  of  this 
effort. 

3.2  Model  Calculations 

In  Figure  3.1  we  compare  electron  profiles,  n(x),  for  several  different 
GaAs  substrates  described  in  Table  3.1.  To  display  electron  profiles  in  a  way 
that  is  comparable  to  the  original  donor  profiles,  Nd(x),  Poisson. s  equation  is 

solved  [25]  using  the  boundary  condition 


^^  =  0atx-0  V 
dx  ’^max* 


(3.1) 


where  x^^ax  is  the  maximum  depth  (2  pm  in  Fig.  3.1)  used  in  the  calculation.  In 
reality,  with  a  Schottky  barrier  at  the  surface  having  a  typical  intrinsic 
barrier  height  of  0.78  V,  n(x)  should  vanish  at  the  surface  when  no  bias  or 
when  reverse  bias  is  applied.  The  boundary  condition,  Eq.(3.1),  allows  defering 
the  choice  of  biasing  conditions  to  a  later  time,  and  reflects  the  physical 
property  that  electrons  do  not  leave  the  surface  nor  the  interior  boundary.  At 
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Substrate  Impurity  Effects  on  Electron  Profiles 


Fig.  3.1.  Calculated  substrate-impuriw  effects  on  electron  profiles 

in  a  channel  of  MESFETs  implanted  using  5x10^2  gj 

ions/cm^.  The  various  substrates  modeled  nave  assumed 

depth-independent  impurity  concentrations  listed  in  Table 

small  depths,  the  electron  concentrations  are  a  factor  of  Tj  =  0  .  8 

below  the  donor  concentrations,  except  at  the  surface,  where  Tl  is  the  assumed 
activation.  The  electron  profile  near  the  surface  is  not  meaningful,  as  it  will 
usually  vanish  under  normal  bias  conditions.  The  donor  profile  for  100  keV  Si 
ions,  used  in  Fig.  3.1,  was  measured  using  secondary-ion  mass  spectroscopy 
(SIMS).  Although  the  wafer  was  oriented  to  avoid  implanting  into  channels, 
the  donor  profile  is  not  Gaussian,  but  shows  an  exponential  tail,  due  to  ions 

scattered  into  channels  in  any  of  the  hundreds  of  collisions  the  Si  ion 

undergoes  before  it  stops.  The  ion  dose  used,  5x10^ ^  cm'^,  was  chosen  to  be 
typical  for  a  microwave  MESFET. 

The  results  shown  in  Fig.  3.1  can  be  understood  in  terms  of  two  factors: 
the  magnitude  of  the  background  acceptor  concentration  and  the  substrate 

barrier  height.  This  is  illustrated  in  Fig.  3.2  where  we  show  dopant 

concentrations  and  band  diagrams  corresponding  to  cases  where  the  substrate 
Fermi  level  is  pinned  at  mid-gap  by  excess  EL2  or  Cr,  and  for  cases  where  the 
Fermi  level  is  pinned  by  excess  shallow  acceptors.  Near  the  intersection  of  the 
implanted  shallow  donor  and  background  shallow  acceptor  in  Fig.  3.2a,  a  p-n 

junction,  called  the  substrate  barrier,  is  formed.  In  the  abrupt  depiction 


Fig.  3.2.  (a)  Dopant  concentrations  and  abrupt-depletion  electron 

profiles  vs  depth.  The  corresponding  energy-band 
diagrams  are  shown  for  cases  where  the  deep  donor  or 
acceptor  concentrations  exceed  the  shallow  dopant 
concentrations  in  the  substrates  (c),  and  for  cases  where 
the  substrate  is  shallow-p-typc  below  the  channel  (b). 

model,  n(x)  would  follow  the  active  donor  concentration  up  to  the  shallow  side 
of  the  substrate  barrier,  then  it  would  fall  off  vertically.  The  Poisson  solutions 
give  approximately  this  behavior,  except  that  n(x)  does  not  decrease  abruptly, 
falling  instead  with  a  typical  Debye  tail  into  the  substrate  depletion  region.  As 
the  background  acceptor  concentration  increases,  the  intersection  between 
the  implanted  donor  and  acceptor  profiles  moves  closer  to  the  surface, 
therefore  the  electron  concentration  profile  is  shallower.  The  other  factor  is 
the  substrate  barrier  height,  which  changes  from  about  0.7  V,  when  the 
substrate  Fermi  level  is  pinned  at  mid-gap  by  EL2  or  Cr  to  about  1.2  V,  when  the 
substrate  Fermi  level  is  pinned  by  excess  shallow  acceptors.  In  the  latter  case, 
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the  large-  substrate  barrier  height  increases  the  width  of  the  substrate 
depletion  region,  so  that  the  electron  profiles  for  Vsub~1.2  V  are  narrower 
than  those  for  V.  In  our  calculations,  we  assume  that  the 

representative  dominant  shallow  acceptor  is  C,  shallow  donor  is  Si,  and  deep 
acceptor  Cr.  Other  acceptor  concentrations  of  Fe,  Mn,  and  Cu,  having  different 
energy  levels  may  be  present  in  LEG  GaAs  [17,19],  which  can  give  other 
barrier  heights  between  0.7  and  1.2  V. 

Consider  the  electron  profiles  for  a  typical  LEG  substrate  where 
the  Fermi  level  is  pinned  by  EL2.  If  we  assume  that  deep  acceptors  are  absent, 
the  shallow'  donor  concentration  is  small,  the  shallow  acceptor  concentration 
is  3.xl0^-  cm -3  and  the  EL2  concentration  is  2x10*^  cm'^  (Table  3.1).  The 


Table  3.1.  Modeled  Substrates  and  Galculated  Threshold  Voltages. 


Substrate 

Na 

Nd 

EL2 

Cr 

Vth 

AVth 

Concentrations  in 

10^®  cm"® 

(V) 

(mV) 

1.  EL2  pinned 

0.3 

0.05 

2 

0 

-6.26 

320 

2.  Heavy  EL2 

03 

0.05 

5 

0 

-6.27 

600 

3.  C  pinned 

0.3 

0.05 

0 

0 

-6.11 

- 

4.  Heavy  C 

1 

0.05 

0 

0 

-5.29 

- 

5.  Cr  pinned 

0.1 

0.2 

0 

1 

-5.59 

190 

6.  Heavy  Cr 

0.1 

0.2  - 

0 

5 

-3.53 

300 

7.  Cr  in  LEC 

0.3 

0.05 

2 

1 

-5.32 

340 

8.  HP  LEC 

0.03 

0 

2 

0 

-6.77 

350 

intersection  between  the  shallow  donor  and  acceptor  profiles  occurs  at  0.38 
pm,  and  the  donor  side  of  the  substrate  depiction  region  extends  about  0.1pm  to 
shallower  depths.  Therefore,  for  x>0.25  pm,  n(x)  falls  off  exponentially. 
Increasing  the  magnitude  of  the  EL2  concentration  has  no  effect  on  the 
electron  concentrations.  Most  of  the  EL2  is  un-ionized,  except  for  the  part 
equal  to  the  shallow  acceptor  concentration,  Na  =  3x10*5  cm'3.  The  Fermi  level 
and  therefore  the  substrate  barrier  height  is  fixed  at  about  0.7  V,  independent 
of  the  EL2  concentration  for  any  concentration  much  above  the  shallow 


acceptor.  This  suggests  that  the  EL2  concentration  should  have  only  a  minor 
effect  on  MESFET  device  characteristics,  unless  the  concentration  falls  below 
the  background  acceptor  level,  which  could  be  caused  by  outdiffusion  [21],  For 
this  case,  the  electron  profile  would  follow  that  for  the  C-pinned  substrate  case 
in  Fig.  3.1.  There,  the  substrate  barrier  height  is  1.2  V,  and  n(x)  falls  off  at 
.shallower  depths  than  where  the  substrate  is  EL2  pinned.  Increasing  the  C 
concentration,  causes  the  point  of  intersection  between  the  shallow  donor  and 
acceptor  to  be  closer  to  the  surface,  resulting  in  a  narrower  electron  profile. 
This  is  found  for  the  high  carbon  concentration  case  in  Fig.  3.1,  where  the 
shallow  acceptor  concentration  is  assumed  to  be  1x10^®  cm'^,  and  the  EL2,  Cr, 
and  shallow  donor  concentrations  are  assumed  to  be  either  zero  or  negligible. 

For  the  deep-acceptor-pinned  substrate,  if  we  assume  that  the  Cr 
concentration  is  1x10^^  cm'3,  and  the  background  shallow  donor  exceeds  the 
background  shallow  acceptor,  as  in  typical  HB-grown,  Cr-doped  GaAs  the 
substrate  barrier  is  formed  between  the  shallow  donor  and  deep  acceptor.  The 
higher  Cr  concentration  used  in  the  Cr-pinned  case  (1x10^^)  shifts  the 
substrate  junction  to  shallower  depths,  resulting  in  n(x)  falling  off  faster. 
Also,  the  electron  concentrations  are  reduced  by  the  background  acceptors. 
This  is  best  seen  in  the  case  of  a  high  concentration  Cr  substrate,  where,  even 
at  the  peak  of  the  n-doping  one  obtains  significantly  reduced  electron 
concentrations  by  acceptors.  For  high-Cr  doping,  the  p-n  junction  is  formed 
at  the  most  shallow  depth,  so  that  the  narrowest  electron  profile  is  obtained 
explaining  the  improved  performance  in  microwave  FETs  fabricated  on  Cr- 
doped  substrates.  Finally,  for  Cr-doped,  LEC-substrates,  n(x)  falls  off  just 
slightly  faster  than  the  case  of  Cr-pinned  substrates.  The  effect  of  the  EL2 
present  in  the  Cr-in-LEC  substrates  is  unimportant,  since  the  Fermi  level  is 
pinned  at  mid-gap  by  the  Cr,  so  EL2  does  not  play  a  role.  A  slightly  higher 
total  acceptor  concentration  is  assumed  for  the  LEC  substrates,  than  for  the  Cr 
substrates.  Note  that  the  electron  profiles  are  steeper  for  the  high-C 
substrates  (total  acceptor=10^ ^  cm’^)  than  for  the  Cr-pinned  ones  (total 
acccptor=l.lx due  to  the  larger  substrate  barrier  height  for  the  C-pinned 
substrates.  This  suggests  that  for  microwave-device  applications,  where  steep 
electron  profiles  are  required,  it  could  be  more  advantageous  to  rely  on  buried 
p  layers  to  sharpen  the  channel  electron  profile  than  Cr  doping  [23].  We  have 
not  considered  Cr  outdiffusion  in  Fig.  3.1  [20],  which  is  difficult  to  control,  and 
can  also  lead  to  changes  in  threshtold  voltages.  For  MMIC  applications, 
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there  should  be  advantages  in  using  the  easily  controlled  buried  p  layer 
instead  of  Cr  [23]. 

Finally  we  consider  the  case  of  implantation  into  high-purity  (HP)  LEC 
substrates.  In  this  case  the  p-n  junction  moves  to  very  large  depths,  so  that 
n(x)  follows  Nd(x)  to  very  large  depths.  We  assume  that  the  background 
acceptor  level  will  still  exceed  the  background  donor  level,  though  both  are 
small.  If  the  background  shallow  acceptor  did  not  the  exceed  shallow  donor, 
the  deep-donor  EL2  could  not  compensate  the  background  impurity,  and  one 
could  not  obtain  a  semi-insulating  substrate.  For  microwave  device 
applications,  slowly  falling  electron  profiles  like  those  obtained  for  HPLEC 

substrates  give  poor  transconductances  at  low  gate  biases  [24].  We  therefore 

conclude  that  background  acceptor  impurity  concentrations  greater  than 
10 cm‘3  are  desirable,  but  they  must  be  uniform. 

Table  3.1  also  gives  calculated  threshold  voltages  for  MESFETs  implanted 
with  5x10^^  100-keV  Si  ions/cm^.  Roughly  speaking,  the  pinchoff  voltage,  Vp> 
is  proportional  to  the  average  depth  of  the  electron  distribution,  so  that  the 
threshold  voltage  Vth  (-Vbi-Vp,  Vbi~0.78  V)  should  be  more  negative  for  deeper 
electron  distributions.  By  comparing  the  electron  profiles  in  Fig.  3.1  and  the 
threshold  voltages  in  Table  3.1,  we  see  that  this  is  indeed  the  case,  Vfh  varies 
from  -3.53  to  -6.77  V  as  one  goes  from  the  steepest  electron  profile  (Heavy  Cr) 
to  the  case  where  n(x)~Nd(x)  (HP  LEC). 

At  small  depths  for  negative  gate  biases  near  threshold,  the  deep  levels 
arc  pulled  above  the  Fermi  level.  The  neutral  deep  donors  are  ionized, 
releasing  electrons  into  the  channel,  and  the  deep  acceptors  neutralize, 

releasing  electrons  into  the  channel.  In  either  case,  the  voltage  needed  to 
pinch  off  the  channel  increases.  The  calculated  relative  shifts  in  the 
threshold  voltages  listed  in  Table  3.1  are  generally  less  than  15%,  even  for  the 
highest  EL2  and  Cr  concentrations.  Such  shifts  are  potentially  observable  in 
threshold-voltage  temperature  dependences,  due  to  the  strong  dependence  of 
the  rale  of  trap  emptying  on  temperature  [25].  The  calculated  shifts  assume 
that  that  the  deep  level  concentrations  are  independent  of  depth.  Due  to  the 
outdiffusion  of  both  Cr  [20]  and  EL2  [6],  the  magnitudes  of  the  deep  level 
concentrations  at  shallow  depths  may  be  very  small,  so  that  much  smaller 
threshold-voltage  shifts  should  be  observed. 


Figure  3.3  shows  the  calculated  effect  of  varying  C  concentrations 
on  MESFET  threshold  voltages.  Carbon  is  usually  measured  by  local- 
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Fig.  3.3.  The  calculated  dependence  of  MESFET  threshold  voltages  on 

the  background  snallow  acceptor  (C)  or  deep  acceptor  (Cr) 
background  levels.  The  ion  dose  and  activation  was 
adjusted  to  give  a  threshold  voltage  of  -0.5  V  for  a  C  level  of 
3x10*5  and  Cr=0.  As  one  increases  either  the  C  or  Cr 

concentrations,  V^h  becomes  more  positive,  with 
approximately  the  same  slope.  The  buried  p  layer  reduces 
the  dependence  on  the  C  or  Cr  levels. 

vibrational-mode  analysis  [18].  In  recent  years  the  usually  accepted  LVM 
absorption  coefficient  [25]  has  inflated  by  a  factor  of  3  to  4  so  that  the 
commonly  observed  range  of  carbon  concentrations  in  EEC  GaAs  has  changed 
from  0.5  to  4x10'^  [16]  to  0.1  to  1x10*^  cm'5  [ig].  This  background  acceptor  has 
two  effects  on  the  MESFET  threshold  voltages.  In  the  abrupt  depletion 
approximation,  Vp  varies  roughly  as 

f*t 

Vp  ~  xdx  [>?Nd(x)-Na(x)],  (3.2) 

Jo 

where  xt  is  the  position  of  abrupt  depletion  edge.  For  higher  C  concentrations, 
the  substrate  p-n  Junction  moves  to  shallower  depths,  so  that  Xt  and  Vp 
decrease.  Second,  the  background  acceptor  concentration  directly  subtracts 


from  the  active  donor  concentration  to  give  a  smaller  pinchoff  voltage,  and 
therefore  a  more  positive  Vth  value. 

For  carbon  concentrations  varying  from  1  to  6x10*^  cm'^,  significant 
changes  in  the  threshold  voltages  are  obtained,  which  must  be  compensated 
for  from  wafer  to  wafer  to  give  uniform  threshold  voltages. 

Deep  and  shallow  acceptors  have  nearly  equal  effects  on  threshold 
voltages  (Fig.  3.3).  Although  models  of  substrate  compensation  and  resistivity 
in  LEC  GaAs  usually  neglect  deep  acceptors  [18),  Cr,  Fe,  Cu,  and  Mn  acceptors 
are  probably  present,  with  concentrations  ranging  from  10*^  to  10^^  cm'^. 
Though  close  to  the  limits  of  SIMS  sensitivity,  Hobgood  et  al.  [17]  reported 
levels  of  Mn  (<1.2x10^^),  Cr  (<1x10'^),  and  Fe  (<mid  10^^  cm‘^)  for  LEC  GaAs 
pulled  from  pyrolitic  Boron  nitride  crucibles.  Kirkpatrick  et.al.  [19]  report 
similiar  levels  (Mn~lxl0^^,  Cr<v5xl0^^,  Fe~3xl0^5).  If  these  concentrations  also 
vary  from  wafer  to  wafer,  and  are  all  electrically  active,  significant  effects  on 
MESFET  threshold  voltages  such  as  those  shown  in  Fig.  3.3  should  be  observed. 

3.3  Conclusions 

Background  impurities  of  shallow  and  deep  acceptors  in  GaAs  substrates 
play  an  important  role  in  both  determining  the  shape  cf  the  electron  profiles 
in  MESFET  channels  and  in  giving  improved  long-term  threshold-voltage 
uniformities  in  the  GaAs  MESFET  industry.  Wafer  manufacturers  should  strive 
to  reduce  impurity  levels  as  much  as  possible.  Uniformity  of  the  impurity 
mu.st  be  achieved  to  increase  MESFET  parameter  control  on  a  wafer  to  wafer 
and  boule  to  boulc  basis. 


IV.  ION-IMPLANTATION  INTO  GaAs  PROFILE  DETERMINATION 


4.1  Ion  Implantation  Profile  Determination  in  GaAs 
4.1.1  Background 

For  controlled  profiles  when  ion  implantating  into  crystalline 
materials,  the  wafers  must  be  tilted  and  rotated  precisely  with  respect  to  the 
beam  to  minimize  axial  and  planar  channeling  (26-33).  This,  however,  does 
not  completely  eliminate  channeling  effects,  since  the  ions  undergo  hundreds 
of  collisions  before  they  come  to  rest,  and  have  a  finite  probability  of 
scattering  into  a  channel  during  any  single  collision.  In  a  channel,  the  ion 
stopping  power  is  reduced,  resulting  in  a  fraction  of  the  ions  coming  to  rest  at 
depths  deeper  than  would  occur  in  an  amorphous  material.  This  partial 
channeling  is  responsible  for  the  exponential  tails  observed  in  implanted 
atomic-impurity  profiles  in  crystalline  materials. 

Several  methods  are  available  to  calculate  profiles  of  ions  implanted 
into  amorphous  materials  (34-38).  Figure  4.1  compares  a  measured  profile 


Fig.  4.1.  Comparison  of  a  measured  80-keV  Si^^  implant 
concentration  profile  with  a  profile  calculated  using  the 
TRIM  Monte-Carlo  code  (35),  and  with  a  fit. 


for  80-keV  Si  ions  implanted  into  GaAs  with  a  profile  calculated  using  the  TRIM 
Monte-Carlo  code  (35)  for  an  identical  implant  into  amorphous  GaAs.  Similar 
results  are  obtained  when  comparing  with  other  implant  calculations;  (34,  36- 
38)  the  amorphous  codes  predict  a  narrower  distribution  and  fail  to  predict  the 
exponentially  falling  channeling  tails. 

In  this  effort  development  of  techniques  to  simulate  ion  implantation 
profiles  in  GaAs,  was  pursued.  GaAs  metal-semiconductor  field-effect 
transistor  (MESFET)  device  chaiacteristics  depend  critically  on  the  threshold 
voltage,  Vth’  which  in  turn  depends  on  the  pinch-off  voltage,  Vp>  given  by 

pOO 

^  xdx  Nd(x), 
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(30)  where  Vbi  (-0.8  V)  is  the  barrier  height  of  the  Schottky  gate  junction. 
The  pinch-off  voltage  is  directly  proportional  to  the  first  moment  of  the 
implanted  donor  concentration  distribution  Nd(x).  As  a  result  of  channeling, 
a  fraction  of  the  ions  are  implanted  at  deeper  depths,  resulting  in  higher 

moments,  and  pinch-off  voltages,  and  more  negative  threshold  voltages,  then 
calculated  for  equal-dose  implants  into  amorphous  material.  Since  quantities 

like  the  drain  source  current  at  zero  gate  voltage  depend  upon  the  square  of 

the  threshold  voltage,  it  is  vital  to  correctly  model  implant  profiles  into 
crystalline  GaAs  to  be  able  to  correctly  predict  device  characteristics. 

The  method  we  adopted  to  simulate  ion-implantation  profiles  in  GaAs  is 
to  fit  analytical  distributions  to  measured  profiles.  The  advantage  of  this 
method  is  that  numerical  profiles  can  be  calculated  rapidly,  as  opposed  to 
Monte-Carlo  or  Boltzmann-transport  calculations  (35-36).  Also,  presently, 

Monte-Carlo  and  Boltzmann  transport  codes  incorporating  crystalline 
structures  have  not  been  sufficiently  developed  to  be  widely  available  (40-44). 
The  disadvantage  of  our  approach  is  that  one  can  never  fully  explore  all 
implant  profiles  of  conceivable  use.  In  this  work,  we  discuss  profiles  for  Be, 
Si,  and  Se  ions  in  GaAs.  These  ions  were  chosen  because  Si  and  Se  are  the  most 
widely  used  donor  impurities  in  GaAs  MESFET  technology,  and  Be  is  used  to 
produce  buried  acceptor  layers  in  MESFETs.  Furthermore,  these  three  ions 
encompass  a  wide  range  of  possible  atomic  numbers,  so  that  eventually  it  may 


be  possible  to  interpolate  parameters  for  calculating  the  profiles  of 
intermediate  atomic  numbers  using  data  from  these  three. 

In  Sect.  4.1.2  we  describe  the  secondary-ion-mass-spectrometry  (SIMS) 
measurements.  In  Sect.  4.1.3,  results  of  profile  studies  are  described  for  (A)  the 
ion-energy  dependence  in  bare,  uncapped  GaAs  wafers,  (B)  the  wafer  tilt  and 
rotation  dependence,  (C)  the  dependence  on  the  ion  dose,  and  (D)  on  the  wafer 
dislocation  density.  In  Sect.  4.1. 3. E  methods  are  formulated  for  calculating  the 
profiles  of  ions  implanted  through  Si02  or  Si3N4  caps,  and  in  Sect.  4.1,3.F 
measurements  of  the  profiles  for  "knock-on"  Si  and  0  atoms,  resulting  from 
implantation  through  Si02  caps  are  discussed.  Section  4.1.3.G  comments  on  the 
relationship  between  SIMS  measurements  of  donor  profiles  and  electron 
carrier  profiles  in  liquid-encapsulated-Czochralski  (LEG)  grown  GaAs.  Finally, 
Sect. 4. 1.4  contains  the  conclusions. 

4.1.2  Experimental  Measurements 

Samples  of  commercial,  high-pressure,  LEC-grown,  (lOO)-oriented  GaAs 
wafers  were  implanted  with  Be^,  Si^^,  and  Se^®  ions  at  doses  of  approximately 
10^3  cm '2  (except  where  noted).  SIMS  measurements  of  the  as-implanted 
profiles  were  made  at  Charles  Evans  and  Associates.  This  ion  dose  was  chosen 
to  obtain  adequate  SIMS  statistics  for  modeling.  Unless  specified,  the  wafer  tilt 
and  rotation  angles  were  7°  and  9°  from  the  wafer  major  flat,  respectively. 
The  samples  were  cleaved  parallel  to  the  major  (Oil)  flat,  and  were  bonded  to  Si 
wafers  at  the  desired  rotation  for  implantation,  or  were  implanted  using  a 
goniometer  to  position  the  samples. 

For  each  measured  SIMS  profile,  a  constant  background  concentration, 
obtained  by  averagihg  the  concentrations  in  the  deep,  flat  portions  of  the 
spectra  was  subtracted.  Also,  since  these  profiles  often  show  a  high, 

unphysical  concentration  in  the  first  few  channels,  we  replaced  the  near¬ 
surface  concentrations  Cj  by  values  obtained  by  extrapolating  the  first  few 

realistic  points  in  each  spectrum  to  shallower  depths.  Using  these  corrected 

profiles,  we  then  calculated  the  moments  of  the  distribution  (x")  using 

<4.2, 

‘  i 
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For  Pearson-lype  distributions,  the  normalized  moments  are  defined  as  (45-46) 


and 


Rp  =  (x‘) 

ARp  =  >/(x^)-(x)^' 


(x-Rp)^ 

ARp3 


0  = 


(x-Rp>^ 

ARp^ 


(4.3) 


where  the  first  moment  Rp  is  referred  to  as  the  projected  range,  and  is  the 
position  of  the  peak  of  the  distribution  for  Gaussian  distributions.ARp  is  the 

width  of  the  distribution,  projected  on  the  depth  axis.  The  normalized  third 
moment  y  describes  the  skew  of  the  distribution.  Positive  values  of  the  order  of 

unity  and  larger  describe  distributions  falling  off  approximately 
exponentially  at  large  depths,  while  negative  values  describe  distributions 
skewed  toward  shallower  depths.  The  kurtosis  factor  p  mainly  decides  the  type 

of  mathematical  function  best  describing  the  distribution  according  to 
whether  it  falls  between  certain  ranges  that  are  functions  of  y  (sec  e.g.  Eq.  4.4) 

below)  (45-47). 

We,  are  mainly  interested  in  the  shapes  of  the  distributions.  In  most  of 
the  data  shown  below,  the  profiles  have  been  normalized  to  arbitrary  peak 
heights;  the  absolute  peak  heights  can  be  obtained  by  normalizing  the 

integrated  area  to  the  implanted  ion  doses. 

To  reduce  SIMS  artifacts  and  uncertainties,  we  measured  most  atom 
distributions  at  least  twice,  using  two  different  SIMS  setups.  Also,  some 

distributions  were  independently  implanted  at  several  different  times  using 

up  to  two  different  implanters.  When  differences  were  observed,  different 
samples  were  measured  two  additional  times  using  the  same  SIMS  setup.  These 
precautions  were  taken  to  increase  our  confidence  that  these  data,  especially 
for  the  energy  dependence  of  the  unencapsulated  implants,  are  accurate,  and 
are  independent  of  cither  the  SIMS  measurement  or  the  ion-implant 

apparatus. 


4.1.3  Results 

4. 1.3.1  Ion-Energy  Dependence 
Figures  4.2  and  4.3  show  moments  for  the  distributions  of  Be,  Si, 
and  Se  ions  implanted  into  GaAs  as  a  function  of  the  ion  energy,  while  Figs.  4.4, 
4.S,  and  4.6  compare  the  measured  distributions  with  profiles 

calculated  using  Pearson-lV  distributions.  To  use  a  Pearson-lV 


ION  ENERGY  E  (keV) 

Fig.  4.2.  Measured  moments  Rp  and  ARp  of  concentration 

distributions  of  20  to  400'keV  Be,  Si,  and  Se  ions  implanted 
into  GaAs,  and  power-law  fits  (for  ion  energies  in  keV). 
The  lO^^-ions/cm^  implants  were  done  into  bare  wafers  at 
a  tilt  and  rotation  angle  of  7°  and  9°,  respectively. 


distribution  for  any  set  of  moments,  the  value  of  the  normalized  fourth 
moment  B  must  be  greater  than  Bq,  defined  as  (45-47) 

o  ^  397^^-48-i-6f-y»i-4)3/2  ,4  4. 

■°o  32_.^2 

We  found  that  the  measured  moments  for  the  profiles  typically  satisfy 
this  criterion.  Using  Pearson-IV  distributions  to  predict  the  profiles  of  ions 
implanted  into  crystals  has  the  advantage  of  built-in  exponential  tails.  An 
alternative  way  often  used  for  Fitting  measured  implant  profiles  is  to  add  an 
exponential  tail  to  a  Gaussian  distribution  (46).  which  is  calculated  using  an 
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Fig.  4.3.  Normalized  third  and  fourth  moments  y  and  6  of 

concentration  distributions  of  Be,  Si,  and  Se  ions  implanted 
into  GaAs,  as  in  Fig.  4.2.  Bq  is  defined  in  Eq.  (4.4). 

The  Y  values  were  fit  to  y  =2.38—1.1  InE,  for  Be,  y  = 

minimum  (2.1,  95E'0-93)^  for  gi,  and  y  =  minimum  (2.6,  6.2- 
1.8  InE),  for  Se,  where  E  is  the  ion  energy  in  keV. 

amorphous  implant  code.  Such  distributions  require  the  same  number  of 
fitting  parameters,  but  the  parameters  are  not  as  straightforwardly  obtained 
from  the  distribution  moments  as  they  are  for  Pearson-lV  distributions.  On  the 
other  hand,  the  Pearson-lV  distribution  is  more  complicated  to  compute  than  a 
Gaussian,  and  can  easily  give  computer  overflow  errors  if  not  properly 
calculated  (48). 

The  moments  Rp  and  ARp  shown  in  Figs.  4.1.2  and  4.3  were  fit  to  power- 
law  functions  of  the  ion  energy,  y  was  fit  to  either  a  linear  function  of  log(E) 

or  a  power-law  function  of  E,  with  a  maximum,  and  B  was  fit  io  a  constant 

multiplied  by  Bq.  Due  to  statistical  uncertainties,  the  third  and  fourth  moments 
are  less  accurate  than  the  first  two;  however,  they  also  have 
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Fig.  4.4.  SIMS  measured  profiles  for  20-200keV  Be  ions  in  GaAs, 

compared  with  Pearson-IV  fits  calculated  with  the 
parameters  given  in  Figs.  4.2  and  4.3.  The 
concentration  distributions  are  normalized  to  arbitray 
peak  heights. 


Fig.  4.5.  SIMS  measured  profiles  for  25-400  keV  Si  ions  in  GaAs,  compared 
with  Pearson-IV  fits. 


less  effect  on  the  derived  distributions.  We  have  little  data  to  support  the  use 
of  a  constant  y  value  for  Si  and  Sc  ions  below  100  kcV;  the  present  form  was 
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Fig.  4.6.  SIMS  measured  profiles  for  50-200-keV  Se  ions  in  GaAs, 

compared  with  Pearson-IV  fits. 

chosen  for  simplicity,  instead  of  some  more  complicated  parabolic  function  of 
the  log(E). 

4. 1.3.2  Tilt  and  Rotation  Dependence 
We  investigated  the  effect  of  varying  the  wafer  tilt  and  rotation  angles 
Stilt  and  Srot  lo  reduce  the  magnitudes  of  the  channeling  tails.  Ziegler  and 

Levers  (28)  have  mapped  the  channel  directions  in  Si  (essentially  the  same 
lattice  structure  as  GaAs)  using  He-backscattering.  As  a  function  of  the 
rotation  angle,  planar  channeling  should  occur  at  0°,  18°,  and  45°  from  the 
(001)  pole,  and  axial  channeling  should  occur  at  various  tilt  angles  for  a 
rotation  angle  of  22°.  The  map  suggests  that  the  best  rotation  angle,  with  a 
standard  tilt  of  about  7°,  should  be  at  about  9°  or  36°  from  the  major  flat  on 
(Oil)  for  SEMI-standard  (lOO)-oriented  LEG  GaAs.  On  the  other  hand,  the  widest 
planar  channel  is  the  (Oil).  To  avoid  this  channel,  the  impinging  angle  0iio. 

given  by  (31), 

‘=^"^110  =  (4.5) 
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should  be  greater  than  about  twice  the  critical  angle  Yc»  given  by  (31) 


where  Zi  and  Z2  are  the  projectile  and  target  (-'32  for  GaAs)  atomic  numbers,  N 
is  the  target-atom  density,  dp  the  distance  between  planes,  and  a  is  the  Firsov 
screening  length.  For  60  keV  Si  ions,  is  equal  to  ~7°.  For  these  ions  ©no  is 
less  than  twice  the  critical  angle  for  every  tilt  angles  less  than  10®  and  every 
rotation  angle,  suggesting  that  avoiding  planar  channeling  is  impossible  for 
this  implant. 

However,  the  magnitude  of  the  channeling  tails  possibly  can  be  reduced 
by  keeping  the  tilt  and  rotation  angles  as  large  as  possible.  Several  authors 

have  suggested  using  tilt  angles  of  at  least  7®  and  rotation  angles  of  45®  for 

implantation  into  Si  and  GaAs  (26-27,  30-31).  Even  if  0iio  is  larger  than  2\|fc 

this  will  not  completely  avoid  channeling,  since  partial  channeling  results 

from  the  multiple  collisions  the  ion  undergoes  in  the  crystal.  As  an  example, 

for  60-keV  Si  ions,  the  TRIM  Monte-Carlo  code  (35)  calculates  that  after 
penetrating  10  nm  of  GaAs,  the  distribution  of  ion  angles  is  Gaussian  [as  in 

most  theories  of  multiple  scatterings  (49)]  with  a  width  ae=20®.  Therefore, 

although  one  may  use  angles  that  avoid  directly  impinging  on  a  channel, 

scattering  into  any  number  of  axial  or  planar  channels  is  possible  after  the 

ion  enters  the  crystal.  For  high  index  planar  and  axial  channels,  the  ion 
lifetimes  in  the  channel  are  short.  Nevertheless,  this  type  of  channeling  will 

contribute  to  the  channeling  tails.  Therefore  independent  of  the  initial  tilt 

and  rotation  angles  chosen  for  the  implant,  there  will  doubtless  be  some 

minimum  channeling  tail. 

We  made  a  series  of  measurements  at  different  tilt  and  rotation  angles. 
A  simple  goniometer  was  constructed  for  our  implanter,  that  allowed  samples, 
cleaved  parallel  to  the  wafer  flat,  to  be  implanted  at  precisely  selected  tilt  and 
rotation  angles.  The  uncertainties  in  these  angles  are  approximately  ±1®. 
Samples  were  also  implanted  in  a  standard  implantor  at  a  tilt  angle  of  7®.  Here 
the  different  rotation  angles  were  obtained  by  bonding  the  samples  to  4-inch 

Si  wafers  at  the  desired  rotation  angles,  allowing  the  simultaneous 
implantation  of  six  different  rotation  angles.  The  relative  uncertainties  in 
these  rotation  angles  are  less  than  I.S®  (resulting  from  the  accuracy  of  the 

bonding  and  the  beam  scan  angle  over  the  limited  part  of  the  Si  wafer  where 
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the  samples  were  placed).  The  absolute  uncertainty  is  about  2.5°.  due  to  the 
placement  of  the  Si  wafer  in  the  implanior. 

Figure  4.7  shows  SIMS  results,  obtained  for  60-lceV  Si  for  selected  tilt 
angles,  with  a  constant  rotation  angle  of  10°.  Narrower  profiles  are  obtained 
at  large  tilt  angles  for  two  reasons:  the  reduction  of  ion  channeling  and  the 
simple  projection  of  the  distributions  onto  the  depth  direction,  according  to 
CO  s9  tilt-  For  the  small  tilt  angles  shown  in  Fig.  4.7,  the 


pjg_  4.7,  Measured  concentration  profiles  for  60-keV  Si  ions 

implanted  into  bare  GaAs  wafers  at  tilt  angles  from  2  to 
11°.  for  a  constant  rotation  angle  of  10°  (except:  0roi=2  6° 
for  9tili=n°)-  The  measured  profiles  are  smoothed  for 
clarity. 

projection  narrows  the  profiles  by  only  0.98,  which  is  clearly  a  negligible 
effect  compared  to  the  observed  narrowing  between  2°  and  11°.  The  broader 
width  of  the  profiles  at  low  angles  is  due  to  a  combination  of  axial  md  planar 
channeling.  One  cannot  completely  avoid  the  channeling  tails;  for  0>1O°,  the 

profiles  are  nearly  independent  of  the  tilt  angle  (aside  from  the  projection 
effect),  but  do  not  approach  the  Gaussian,  predicted  for  implantation  into 
amorphous  material.  Similar  results  were  obtained  for  100-  and  240-keV  Si  ion 
implantation  into  GaAs.  (26,31)  and  for  Implantation  of  various  ions  into  Si 
wafers  (29,  31-33). 
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In  Fig.  4.8  we  plot  the  rotation-angle  dependence  for  20-kcV  Be 
implantation  into  GaAs,  for  a  tilt  angle  of  7°.  Since  in  this  case  2\\f  q  is  ~11®, 

channeling  into  the  (Oil)  plane  cannot  be  avoided  for  any  rotation  angle  at 
this  tilt,  although  slightly  narrower  profiles  can  be  obtained  at  large  rotation 
angles,  as  is  clearly  seen  in  the  figure. 


0.0  0.1  0.2  0.3  0.4 
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Fig.  4.8.  Measured  concentration  profiles  for  20-keV  Be  ions 

implanted  into  GaAs  for  selected  rotation  angles  from  0®  to 
45°,  for  a  constant  tilt  angle  of  7°. 

For  200-keV  Be  ion  implantation,  2\fc  is  3°.  With  0tiii=7°,  ®110  is  greater 
than  2\|rc  for  rotation  angles  greater  than  10°,  significantly  reducing  (011) 
planar  channeling.  As  is  seen  in  Fig.  4.9,  the  profile  for  0°  rotation  is  much 
broader  than  those  obtained  for  all  other  angles,  in  agreement  with  this 
expectation.  In  several  cases,  we  purposely  implanted  at  rotation  angles  to 
enhance  channeling  into  minor  planes  mapped  by  Ziegler  and  Lever  (28). 
Aside  from  the  (011).  the  most  significant  should  be  the  (400)  at  a  rotation 
angle  of  45°.  However,  impinging  on  these  minor  planes  is  seen  to  have  no 
measurable  effect  on  the  channeling  tails;  identical  profiles  are  seen  at 
nearby  rotation  angles. 

By  implanting  through  amorphous  overlayers,  a  reduction  in  the 
dependence  of  the  implant  profiles  on  the  tilt  and  rotation  angles  can  be 
achieved  (26).  The  overlayer  scatters  the  beam  sufficiently  to  reduce  the 
fraction  of  ions  directly  impinging  upon  the  axial  or  planar  channels.  For 
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Fig.  4.9.  Measured  concentration  profiles  for  200-keV  Be  ions 

implanted  into  bare  and  capped  (45-nm  Si2N4)  GaAs  wafers 
for  selected  rotation  angles  from  0°  to  45°  and  at  a  constant 
till  angle  of  1°.  The  wide  black  lines  indicate  the  envelope 
of  the  data  obtained  for  five  rotation  angles  between  9° 
and  45°. 


example,  the  difference  between  the  profiles  resulting  from  the  0rot=O° 
implant  and  other  rotation  angles  for  200  keV  Be  ions  implanted  through  a  45- 

nm  cap  Si3N4  layer  is  smaller  than  that  seen  for  the  bare-wafer  implants  (Fig. 

4.9). 

4. 1.3. 3  Dose  Dependence 

Due  to  the  lattice  displacements  produced  by  the  stopping  of  heavy  ions 
in  crystalline  GaAs,  one  can  expect  to  heavily  damage  the  crystal  in  the  region 

where  the  ion  comes  to  rest.  At  high  doses,  a  profile  that  is  closer  to  Gaussian 

or  Pearson-1  profiles,  predicted  by  the  amorphous  implant  codes  should  then 
be  obtained  (35).  On  the  other  hand,  high  implant  ion  currents  can  cause  the 
wafer  temperatures  to  rise,  resulting  in  some  self  annealing  taking  place,  thus 
repairing  the  damage,  and  resulting  in  the  same  relative  channeling  tails 
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being  obtained  for  both  low  and  high  doses.  In  some  Si  device  fabrication 
processes,  the  crystal  is  purposely  heavily  damaged  to  suppress  channeling. 
Producing  heavy  damage  in  GaAs  is  not  desirable,  since  residual  defects 
usually  remain  following  the  post-implant  anneal,  resulting  in  electrically 
active  traps,  low  ion  activation,  and  poor  electron  mobilities. 

To  be  of  use  for  the  fabrication  of  MESFET  devices,  the  profiles 
calculated  with  our  analytical  formulas  should  be  valid  for  doses  between  10^  ^ 
and  10^3  ions/cm^,  typically  used  for  the  MESFET  active  channel.  Doses  higher 
than  lO^^cra'2  are  sometimes  used  to  implant  the  source  and  drain  contact 
regions  of  MESFETs,  but  device  characteristics  are  not  critically  dependent  on 
the  shape  of  those  profiles.  To  assure  that  the  profiles  are  independent  of  ion 
dose  below  the  10^^  dose  we  used  for  most  of  our  work,  implants  of  3x10^2  jq 
10^5  ions/cm^  were  done.  Figure  4.10  shows  profiles  for  100  keV  Si  ions. 


Fig.  4.10.  Measured  concentration  profiles  for  100-keV  ^^Si  ions 

implanted  into  bare  GaAs  wafers  at  doses  of  10^  3,  lO^'*,  and 
10*^  ions/cm^.  The  profiles  have  been  normalized  to  the 
same  peak  height,  2x10*^  cm'^. 

normalized  to  same  peak  concentration.  At  higher  doses,  narrower  profiles 
are  observed.  Yeo  et.al.  (50)  obtained  similar  results  for  120-keV  Mg  ions  in 
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GaAs,  and  we  obtained  similar  results  for  40  and  200-keV  Be  and  25-,  60-  and 
200-keV  Si  ion  implants.  The  reduction  of  the  channeling  tails  is  smaller  for 
Be  than  for  the  more  heavily  damaging  Si  ions.  Consistent  results  are  obtained 
below  10^^  ions/cm^,  so  that  the  results  given  in  Sect.  4.1  for  profiles 
measured  for  10^^  ions/cm^  doses  should  be  valid  for  implant  doses  of  10*2 
cm  '2. 

Selenium  produces  the  most  damage,  and  Se-ion  implantation  showed 
the  greatest  effect  of  increasing  dose  on  the  profiles,  as  seen  in  Fig.  4.11. 
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Fig.  4.11.  Normalized,  measured  concentration  profiles  for  200-keV 

Se  ions  implanted  into  bare  GaAs  wafers  at  doses  of  10*2, 
10*"^,  and  10*5  ions/cm2.  Fitted  Pearson-IV  profiles,  and  a 
TRIM  (Ref.  35)  calculated  profile  for  implantation  into 
amorphous  GaAs  are  shown  for  comparison. 


Between  10*2  and  5x10*2  ions/cm2  the  magnitudes  of  the  channeling  tails 
decreased,  with  the  reduction  beginning  to  saturate  above  lO*'*  cm '2.  Despite 
the  very  high  doses  and  resultant  damage,  we  never  obtained  profiles  close  to 
those  for  implantation  into  amorphous  GaAs  (35).  This  may  be  due  to  self¬ 
annealing  effects  associated  with  wafer  healing  for  the  high  doses.  (In  the 
high-current  Se  implants,  the  wafer  temperatures  became  sufficiently  high  to 
partially  melt  the  photoresist  used  to  bond  the  samples  onto  Si  wafers.) 
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Figure  4.11  suggests  that  to  obtain  narrower  profiles,  a  room- 
temperature  damage  implant  with  As  (being  similar  to  Se)  at  5x10*^  ion/cm^  is 
sufficient.  The  additional  damage  may  result  in  high  post-implant  trap 
concentrations,  (51)  though  there  is  some  evidence  that  using  As  ions  may 
result  in  increased  activation  for  subsequent  Si  implants  (52). 

4. 1.3. 4  Dislocation  Density  Dependence 

The  commercial,  high-pressure-grown  GaAs  wafers  used  in  our  studies 
have  typical  dislocation  densities  of  3  to  8x10^  cm'^  (53).  In  a  typical  SIMS 
measurement,  the  sputter  crater  diameter  is  about  100  llm,  so  approximately  6 

dislocations  are  expected  to  be  present  in  every  measurement.  To  investigate 

whether  these  dislocations  might  scatter  the  implant-ion  beam  sufficiently  to 
affect  the  channeling  tails  (54),  we  compared  profiles  for  low-dislocation.  In- 
doped  wafers  and  conventional  wafers  (55),  both  implanted  with  10^^  60-keV 
Si,  200-keV  Se,  and  100-keV  Be  ions/cm^.  For  the  (2%)  In-doped  wafers, 
varying  dislocation  densities  are  obtained  in  the  middle  of  the  wafer  (3000 
cm ’2),  in  a  ring  approximately  l-inch  in  diameter  on  a  3-inch  wafer  (500 
cm-2),  and  on  the  outer  part  of  the  wafer  (800  cm*^).  We  also  compared 
conventional  wafers  having  a  dislocation  density  equal  to  45000  cm*^.  As  seen 
in  Fig.  4.12,  identical  implant  profiles  are  obtained  with  all  dislocation 

densities  for  Se.  Similar  results  are  obtained  for  Si  and  Be  implantation. 

We  conclude  that  the  results  obtained  in  Sect. 3.1  for  implantation  into 
high-dislocation  material  should  also  be  valid  for  implantation  into  low- 
dislocation,  In-doped  material. 

4. 1.3. 5  Implantation  Through  Caps 

In  the  fabrication  of  GaAs  devices,  implantation  through  40-  to  100-nm 

Si3N4  or  Si02  encapsulants  is  sometimes  done.  These  caps  prevent  the 
decomposition  of  the  GaAs  during  subsequent  post-implant  annealing,  and  also 
prevent  the  deposition  of  contaminants  onto  the  surface  and  possible  surface 
decomposition  resulting  from  wafer  heating  during  the  implant  and  other 

processing.  We  have  seen  in  Sect.  3.2  that  implanting  through  caps  reduces 
the  dependence  of  the  channeling  tails  on  the  wafer  rotation  angle,  a  variable 
not  often  accurately  controlled  in  many  manufacturing  facilities. 

We  have  measured  the  profiles  for  Si,  Se,  and  Be  ions  in  GaAs,  implanted 
through  40  to  200-nm  plasma-enhanced  chemical-vapor-deposition  (PECVD) 
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200  keV  Se  in  GaAs 
Dislocation  Density  Dependence 


DEPTH  (pm) 

Fig,  4.12.  Measured  concentration  profiles  for  200-keV  Se  ions 

implanted  into  bare  GaAs  wafers  with  selected  dislocation 
densities,  obtained  by  using  In-doped  GaAs  wafers  (etch- 
pit  density  -  500,  800,  and  3000  cm’2)  and  conventional 
wafers  (EPD  =  45000  cm*^). 

Si3N4  and  CVD  Si02  layers.  The  estimated  cap  densities  are  2.9  and  2.2  g/cm^, 
respectively. 

For  Be  ions,  implanting  through  caps  mainly  just  shifts  the  origin  of 
the  implant  profiles.  The  range  scaling  theory  (45-47)  was  developed  to 
describe  the  concentration  distributions  both  in  the  cap  Ccap(x)  and  the 
semiconductor  Cs-  Since  the  caps  are  amorphous,  TRIM  (35)  or  other 
amorphous  profile  simulators  can  be  used  to  calculate  those  concentration 
distributions.  In  the  underlying  semiconductor,  the  concentration  is  given  by 

C(.x)  =  C5(x-tc3pRp5/Rpcap)  (4-. 7) 

where  x  is  measured  from  the  surface  of  the  semiconductor,  and  Rps  and  Rpcap 
are  the  projected  ranges  of  the  ions  in  the  semiconductor  and  in  the 
encapsulant,  respectively.  In  the  present  case,  we  use  the  Pearson-lV  profiles, 
derived  for  implantation  into  unencapsulated  wafers  in  Sect.  3.1,  to  estimate  C^. 
The  amorphous  implants  are  usually  described  analytically  by  Gaussians  (Si 
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and  Se  ions)  or  by  Pearson-I  profiles  (Be  ions).  The  ion  dose  lost  in  the  cap  is 
obtained  by  integrating  over  the  profile  in  the  cap: 


f  *cap 

*cap-  dit' Ccap(x'). 

Jo 

and  the  concentration  in  the  GaAs  is  then  normalized  to  the  total  implant  dose 
*I*minus  *^cap- 

For  the  ion  range  in  the  semiconductor,  Rps»  a  choice  must  be  made 

whether  to  use  the  amorphous  range,  consistent  with  Rpcap*  or  the  larger, 

experimentally  determined  projected  range  shown  in  Fig.  4.12,  which  is 

influenced  by  channeling.  In  fitting  the  measured  profiles  of  Be  ions 

implanted  through  caps,  we  find  that  using  the  amorphous  range  gives  better 
agreement  with  the  data,  as  shown  in  Fig.  4.13.  Similar  results  are  obtained  for 
100-  to  400-keV  Be+  ion  implants  through  50  to  200  nm  Si3N4  and  Si02  caps. 


Fig.  4.13.  Measured  concentration  profiles  for  50-keV  Be  implanted 

into  bare  GaAs  wafers  and  wafers  with  45-  and  79-nm  CVD 
Si3N4  caps,  compared  with  fits  using  the  range-scale 
theory,  Eq.  (4.7). 

Wc  have  also  calculated  the  amorphous  ranges,  and  the  longitudinal  and 
lateral  straggling  parameters  using  the  TRIM  Monte-Carlo  code  (35).  In  this 
calculation,  about  2000  to  5000  ion  trajectories  are  followed,  and  the  ranges  are 
fit  to  power-law  functions  of  the  ion  energy,  e.g.. 
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The  parameters  A  and  B  determined  for  implantation  into  GaAs,  Si3N4 
and  Si02  are  listed  in  Table  4.1.  For  caps  with  different  densities,  one  scales  the 
three  moments  inversely  with  density  .  For  the  densities  listed,  we  generally 
find  that  the  ratio  of  Rp  values  in  Eq.  (4.7)  is  near  unit  for  Si3N4  and  about  2/3 
for  Si02  encapsulants. 

TABLE  4.1.  Amorphous  ranges  Rp,  straggling  ARp,  and  lateral 

straggling  ARi  Power-law  parameter  [Eq.  (35)]^. 


Ion 

Substrate** 

R 

p 

<1 

A/?, 

A 

B 

A 

B 

A 

B 

Si^’ 

GaAs 

13.021 

0.951 

13.461 

0.777 

11.477 

0.831 

Si 

Si,N4 

11.187 

0.983 

9.223 

0.773 

5.371 

0.870 

Si 

SiOj 

13.850 

1.020 

12.349 

0.793 

6.341 

0.919 

Se«° 

GaAs 

6.049 

0.932 

4.698 

0.830 

3.423 

0.860 

Se 

SijN^ 

7.853 

0.902 

3.211 

0.836 

2.319 

0.839 

Se 

SiO, 

11.499 

0.892 

4.893 

0.813 

3.536 

0.822  . 

Mg^" 

GaAs 

15.770 

0.952 

15.898 

0.781 

13.095 

0.851 

Mg 

SijN^ 

13.532 

0.986 

13.175 

0.732 

6.596 

0.881- 

Mg 

SiO, 

17.901 

1.002 

15.381 

0.779 

8.765 

0.903 

Be’ 

GaAs 

54.931 

0.846 

47.144 

0.647 

55.350 

0.694 

Be 

Si3N4 

47.828 

0.882 

39.607 

0.583 

28.388 

0.724 

Be 

SiO, 

69.352 

0.882 

60.177 

0.576 

56.138 

0.585 

*The  moments  are  in  A,  when  the  ion  energy  is  in  keV. 
‘’Densities;  GaAs:  5.32  g/cm^  Si,N4:  3  g/cm^;  SiO,:  2.2  g/cm\ 


The  magnitudes  of  the  channeling  tails  are  reduced  when  Si  and  Se  ions 
are  implanted  through  caps.  This  is  illustrated  in  Fig.  4.14,  together  with  the 

method  used  to  fit  this  reduction.  For  100-keV  Si  implants  through  79  nm  Si3N4 
caps,  simply  shifting  the  Pearson-lV  profile,  obtained  for  implants  into  bare 

wafers,  results  in  a  larger  magnitude  of  the  tail  at  large  depths,  though  it 

accurately  predicts  the  shape  of  the  measured  distributions  at  large 
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Fig,  4.14,  Measured  concentration  profiles  for  100-keV  Si  implanted 

into  GaAs  through  79-nm  CVD  Si3N4  caps  compared  with 
calculated  fits. 

depths.  We  fit  the  measured  distribution  by  reducing  the  magnitude  of  the 
shifted  Pearson-lV  (calculated  using  Eq.  (4.7),  by  a  factor  of  0.25.  For  the 
remaining  part  of  the  distribution,  a  shifted  Gaussian,  calculated  using  the 
TRIM  parameters  in  Table  4.1,  and  Eq,  (4.7)  is  used  to  describe  the  shift.  Both 
the  Pearson-1  V  and  Gaussian  distributions  are  normalized  to  unity  at  x=Rp,  so 
the  mixture  can  be  described  by  a  single  parameter  f 

C(x)=.(  I  -f)CsGaus(x)  +  rcsPIv(*)  (4. 10) 

The  net  concentration  is  then  normalized  to  the  integrated  dose  or  to  an 
arbitrary  peak  height,  as  shown  in  Fig.  4.14. 

Presently,  insufficient  data  exists  to  obtain  a  universal  parameterization 
for  the  parameter  f.  Probably,  f  depends  on  the  ratio  of  the  cap  thickness  to 
Rp  in  the  semiconductor.  Fits  for  other  100-keV  Si  ion  implants  are  shown  in 
Fig.  4.15.  The  Pearson-1  V  fraction  is  seen  to  decrease  as  the  cap  thickness 
increases.  Apparently,  for  the  relatively  deeper  Be  ion  implants,  the  cap 
thickness  is  not  sufficiently  large  compared  to  Rp  to  observe  a  reduction  in  f. 
Also,  electronic  stopping  dominates  the  Be  energy  loss  in  the  caps,  so  the  Be 
ions  are  not  scattered  as  much  as  the  Si  and  Se  ions.  As  more 
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Fig.  4.15.  Measured  concentration  profiles  for  100-keV  Si  implanted 

through  PECVD  Si3N4  and  one  CVD  Si02  cap  into  GaAs.  The 
numbers  give  the  cap  thickness  in  nm,  and  the  numbers 
in  parentheses  give  the  Pearson-IV  mixing  factor  f 
defined  in  Eq.  (4.40). 

data  becomes  available  for  Si  and  Se  implantation,  further  parameterization  of 
this  effect  will  be  done. 

4. 1.3.6  Knock-on  Distributions 
When  one  implants  into  GaAs  through  Si3N4  or  Si02  caps.  Si  and  N  or  0 
atoms  are  knocked  forward  by  the  implant  ion  into  the  GaAs  (36,  56).  The 
"knock-on"  Si  recoils  in  the  GaAs  are  additional  donors,  which  add  to  the 
effective  concentration  of  donors  in  the  GaAs,  while  the  O  recoils  can  either 
form  deep-donor  traps  or  can  compensate  the  implanted  donors. 

It  is  difficult  to  measure  the  knock-on  distributions,  except  for  high 
implant  doses.  In  the  SIMS  analysis  of  GaAs  substrates,  masses  associated  with 
0.  N,  and  Si^*  typically  have  large  backgrounds.  (Mass  28  is  probably  due  to  CO 
molecules.)  We  have  measured  knock-on  distributions  fo.  an  implant  dose  of 
2x10' 5  200-keV  Se  ions/cm^  implanted  through  87-nm  Si02  caps.  This  dose  is 
sufficiently  high  so  we  can  measure  knock-on  Si^^  from  caps  having  the 
normal  4.7%  Si^^  abundance.  The  O  knock-on  distribution  was  measured  by 
taking  the  difference  between  the  O  found  after  implanting  through  a  Si02 
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cap,  and  the  O  background  found  after  implanting  through  a  74-nm  sputtered 
Si3N4  cap.  Similar  Si  knock-on  distributions  were  observed  for  implants 
through  the  Si3N4  cap,  but  N  atom  concentrations  could  not  be  measured. 

Figure  4.16  compares  the  measured  knock-on  distributions  with  TRIM 
calculations.  As  usual,  the  Se-ion  distribution  has  a  large  channeling  tail, 
though  it  is  not  as  large  for  these  high  dose  implants,  as  it  would  be  for  a  low- 
dose  one  (Fig.  4.11).  The  magnitudes  of  the  Si  and  O  knock-on  concentrations 
agree  roughly  with  the  TRIM  calculations  for  shallow  depths,  but  like  the  Se 
ion  implants,  they  exhibit  longer  tails  than  predicted  by  TRIM. 


DEPTH  (pm) 

Fig.  4.16.  Measured  (lines)  and  TRIM  calculations  (points)  of  the 

concentration  profiles  of  Se  ions,  and  recoiled  Si  and  O 
atoms,  after  implantation  of  2x10^5  200-keV  Se  ions/cm^ 
through  an  87-nm  Si02  cap  on  GaAs. 

Since  the  recoil  concentrations  are  largest  at  the  surface  where 
channeling  tails  are  still  not  significant,  we  conclude  that  TRIM  profiles  may 
be  adequate  for  simulating  recoil-atom  profiles  for  MESFET  process  and  device 
modeling. 

4.1.4  Conclusions 

Profiles  measured  for  Be,  Si,  and  Se  implantation  into  unencapsulated 
GaAs  wafers  at  energies  between  20  and  400  keV  have  been  fit  to  Pearson-lV 
distributions,  and  the  resulting  fitting  parameters  were  fit  to  functions  of  the 
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ion  energy  to  derive  simple,  analytical  formulas  for  implant  profiles  for  any 
energy  between  about  10  and  400  kcV. 

Profile  measurements  for  selected  wafer  tilt  and  rotation  angles  show 
that  the  narrowest  profiles  for  a  given  ion  energy  occur  for  tilt  angles  greater 
than  7°,  and  rotation  angles  of  45°.  Implantation  through  wafers  encapsulated 
with  Si3N4  or  Si02  layers  reduces  the  dependence  of  the  profiles  on  the 

rotation  angle. 

The  profiles  are  not  dependent  on  the  wafer  dislocation  density  for 
densities  between  500  and  50,000  cm'^,  nor  on  In  doping  of  the  wafers. 
Channeling  tails  on  the  profiles  can  be  reduced  slightly  for  Si  and  Be  ions  by 
using  high  ion  doses  of  either  the  bombarding  ion  or  a  similar,  but 

electrically-inactive,  damage-producing  ion.  Se  ions  have  the  greatest  effect 
on  reducing  channeling  tails,  though  it  is  unclear  whether  the  possibly 
beneficial  effect  of  reducing  tails  might  be  outweighed  by  undesirable  effects 
of  trap  formation,  low  ion  activation,  and  poor  electron  mobilities. 

For  Be  implants  through  Si02  and  Si3N4  encapsulants  on  GaAs,  the 
profiles  in  the  GaAs  can  be  predicted  by  simply  shifting  the  profile  for 
implantation  into  unencapsulated  wafers  by  an  amount  related  to  the  cap 

thickness.  For  Si  and  Se,  our  measurements  suggest  that  the  channeling  tails 
are  reduced  slightly  by  implanting  through  caps.  We  presently  have 
insufficient  data  to  derive  a  universal  function  describing  this  reduction.  The 
distributions  of  Si  knock-on  ions  after  implantation  through  caps  are 

predicted  accurately  using  TRIM  Monte-Carlo  calculations  at  shallow  depths. 
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4.2  Device  Sensitivity  Resulting  From  Ion  Implantation 

4.2.1  Background 

Ion  implantation  is  the  most  widely  used  process  in  the 
semiconductor  industry  for  introducing  controlled  impurity  concentrations 
into  Si,  GaAs,  and  other  materials.  However,  many  effects  are  associated  with 
ion  implantation  that  require  a  clear,  quantitative  understanding  to  maximize 
the  performance  and  uniformity  of  devices  manufactured  using  technique. 

The  present  work  examines  how  implantation  profiles  can  affect  GaAs 
MESFET  electrical  device  characteristics.  To  illuminate  implantation  effects, 
we  have  used  the  process  and  device  model  discussed  in  Section  V  to  calculate 
device  threshold  voltages,  currents,  and  transconductances.  Our  present 
calculations  mainly  address  the  issues  of  how  to  maximize  the  performance 
and  uniformity  of  GaAs  MESFETs.  In  consideration  of  device  performance,  a 
primary  topic  we  consider  is  how  much  the  shapes  of  the  implant  and 

resulting  electron  carrier  profiles  limit  the  ultimate  transconductances 

obtainable  for  ion  implanted  MESFET's.  In  Sections  4.2.2  and  4.2,2 
transconductance  calculations  for  different  energy  implants  and  MESFET 
biases  are  shown.  Four  possible  effects  on  device  uniformity  related  to 

implanted  channels  are  considered  in  the  remaining  part  of  Sect.  Ill:  the 
effects  of  cap  thickness  and  etch-depth  fluctuations  on  MESFET  threshold- 
voltage  uniformity,  dopant  diffusion  during  post-implant  annealing,  and  the 
effects  of  Si-recoil  ion  distributions  for  implants  through  caps. 

4.2.2  Transconductance 

4.2.2. 1  Profile  Width  Effects 
Narrower  carrier  profiles  give  higher  MESFET  transconductances.  For 
FET's  operating  in  the  saturated-current  regime,  and  fabricated  using  epitaxial 
layers  of  thickness  a,  the  transconductances  vary  as  [64,  Eq.(23b)] 

(4.11) 

where  W  is  the  FET  width,  Vj  is  the  saturated  electron  velocity,  and  Vg  is  the 
gate  voltage.  A  similar  relationship  should  be  approximately  valid  for  ion- 
implanted  MESFET's,  where  a  is  some  effective  implant  layer  thickness  that 
depends  on  the  ion  energy.  To  gauge  the  dependence  of  a  and  therefore  gm  on 
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the  implant  energy,  we  plot  in  Fig.  4.17,  transconductances  for  FETs  with 
constant  threshold  voltage  V),  constant  Vbi,  constant  gate  voltage 

(Vg^OV),  and  width  (W=l  mm).  Although  this  model  uses  saturated  velocities 
that  depend  logarithmically  on  the  absolute  doping,  the  average  saturated 
velocity  is  approximately  1.2xl0^cm/s  for  all  profiles  shown  in  Fig.  4.17., 
therefore  the  transconductance  variations  with  implant  energy  directly 
reflect  the  variations  of  the  effective  layer  thicknesses.  The  threshold 
voltages  are  kept  constant  by  varying  the  ion  doses,  shown  in  Fig.  4.17. 


(a)  Esi  (keV)  (b)  Esi  (keV) 


Fig.  4.17.  Calculated  transconductances  gm  at  zero-gate  voltage  and 

the  dose  needed  to  obtain  a  constant  threshold  voltage  of  - 
0.5  V  for  MESFETs  fabricated  by  implanting  Si  ions  with 
energies  between  10  and  200  keV  into  GaAs.  For  the  solid 
line,  a  Pearson-IV  profile,  based  on  fits  to  measurements 
was  used.  For  the  dashed  and  dash-dotted  lines,  Gaussian 
profiles  were  used,  based  on  TRIM  [35]  and  LSS  [34] 
calculations  of  Rp  and  A  Rp. 

Three  different  types  of  profiles  are  compared  in  Fig.  4.17,  profiles 
derived  from  the  LSS  tables  [34],  TRIM  Monte-Carlo  profiles  [35],  and  measured 
profiles,  fit  to  Pearson-lV  functions  (See  Section  4.1).  The  LSS  tables  predict 
lower  values  of  Rp  than  TRIM,  which  in  turn  predicts  lower  values  than  are 
observed  experimentally.  Since  the  parameter  a  is  roughly  proportional  to  Rp, 
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and  the  pinch-off  voltage  varies  roughly  as  the  product  of  the  dose  and  a,  to 

keep  Vp  and  therefore  Vth  constant,  higher  doses  are  required  when  using  LSS 
profiles  when  compared  to  the  Pearson-lV  profiles  to  obtain  the  same 
threshold  voltage.  With  the  doses  adjusted,  the  resulting  transconductances 
for  the  LSS  profiles  are  larger  than  for  either  the  TRIM  or  measured  profiles. 

None  of  the  transconductances  increase  inversely  with  the  ion  energy 
E.  In  part  this  occurs  because  Rp  only  varies  as  but  other  factors,  such 

as  ARp,  the  profile  shapes,  and  the  absolute  magnitude  of  the  donor  doping, 
compared  to  the  wafer  acceptor  impurity  concentrations  also  play  a  role  in 

determining  the  effective  layer  thickness.  As  the  ion  energy  decreases  from 

200  to  20  keV,  the  transconductances  for  the  measured  profiles  only  increase 

by  a  factor  of  3  to  4.  Between  30  and  10  keV,  g^  only  increases  from  150  to 

225  mS/mm,  much  less  than  a  factor  of  three.  Sugitani  et.al.[59]  noted  a 

similar  slow  increase  in  effective  layer  thicknesses  for  10,  20,  and  30-keV  Si 

implants  in  GaAs,  derived  by  fitting  observed  distributions  to  Gaussian 

functions,  and  relating  a  to  the  first  and  second  moments  of  the  distributions. 

Since  Se  implant  profiles  for  ion  energies  scaled  according  to. 
Ese  =  0-829Esi*-28  are  essentially  identical  to  Si  profiles  implanted  at  an  energy 
Esi.  [57]  the  results  shown  in  Fig.  4.17  may  also  be  applied  to  Se  channel 
implants. 

4. 2. 2. 2  Maximizing  Low-Bias  Transconductances 

For  MESFETs  fabricated  for  small-signal  microwave  amplifier 
applications,  high  transconductances  are  needed  at  low  gate  voltages  (low  DC 
source-drain  currents  Ids)  to  achieve  high  gain  at  low  power  levels.  The  best 
practical  channel  profile  for  achieving  this  is  a  "box"  or  epitaxial-layer 

profile,  where  the  donor-impurity  concentration  falls  off  from  a  constant 
value  to  background  levels  over  a  few  atomic  layers.  Gaussian  implant  profiles 

are  less  desirable,  but  the  observed  implant  profiles,  with  their  exponential 
tails  due  to  ion  channeling,  are  even  less  optimal.  Figure  4.18  compares 
calculated  transconductanccs  as  a  function  of  Ips/^DSS  for  each  of  these 
profiles  (lDSS  =  fDS(Vg=0, Vd=3  V)).  The  box  profile  has  a  constant  doping 
concentration  of  3x10^2  cm"5,  to  a  depth  of  a=0.165  pm.  The  implant  profiles  for 
100-keV  Si  use  a  Gaussian  profile,  derived  from  TRIM  Monte-Carlo  [35] 

calculation,  and  a  Fearson-lV  profile,  derived  from  SIMS  measurements.  For 

both  the  Gaussian  and  Pcarson-lV  profiles,  a  dose  of  SxlO^'^  ions/cm^  is  used, 
with  recess  etching  to  obtain  a  -3  V  threshold  voltage  performed.  The  channel 
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donor  profiles  are  shown  in  Fig.  4.18b.  Figure  4.18a  shows  that  at  a  gate-bias  of 
OIIdSS*  ihe  transconductances  arc  only  70  mS/mm  for  the  Pearson-IV, 


Fig.  4.18.  (a)  Transconductances  of  recess-etched  MESFET's  versus 

the  DC  source-drain  current,  normalized  to  the  zero-gate 
voltage  value.  Calculations  were  made  using  Pearson-lV 
profiles,  Gaussian  profiles,  and  constant-doping  (BOX) 
donor  profiles,  which  are  shown  in  (b). 

but  80  mS/mm  for  the  box  profile.  Due  to  differences  in  the  exact  shapes  and 
magnitudes  of  the  doping  profiles,  the  three  do  not  give  exactly  identical 
transconductanccs  at  high  Id  S  •  If  they  are  normalized  to  the  same 
transconduclance  at  Idss<  the  differences  between  the  transconductances  for 
the  box  and  Pearson-1  V  profiles  become  even  greater. 

Although  the  Pearson-lV  donor  profiles  are  very  broad  when  compared 
to  the  other  two  profiles,  much  of  the  channeling  tail  is  cut  off  by  the 
depletion  region  of  the  substrate  junction  created  by  the  implanted  donors  and 
the  substrate  acceptor  concentration,  assumed  in  Fig.  4.18  to  be  3x10^^  cm"  ^ 
for  LEC  GaAs  wafers.  The  resulting  electron  profiles  are  therefore  more 
Gaussian  than  the  donor  atom  profiles,  which  significantly  improves  the  low- 
bias  transconductances.  Figure  4.19  shows  electron  profiles  and  the 
corresponding  transconductances  for  various  substrate  acceptor 
concentrations  for  devices  implanted  with  5x10^2^  100-keV  Si  ions/cm^,  and 
etched  to  achieve  a  -3  V  threshold  voltage.  For  a  high-purity  wafer  without 


Fig.  4.19.  Transconduciances  of  recess-etched  MESFET’s  versus  the  DC 

source-drain  current  (normalized  to  the  zero-gate-voltage 
value),  calculated  using  simulations  of  measured  implant 

profiles  for  5x10^^  100-keV  Si  ions/cm^,  and  assuming 
various  wafer  acceptor  impurity  concentrations.  The 

corresponding  electron  carrier  profiles  after  recess 

etching  to  obtain  a  -3  V  threshold  voltage  are  shown  in  b. 

acceptor  impurities,  the  electron  and  donor  profiles  are  nearly  identical,  with 

the  transconductance  at  O.IIdSS  being  less  than  60  mS/mm.  Carbon  acceptor 

concentrations  between  1  and  6x10^^  cm*^  can  be  present  in  LEC-grown  GaAs. 
Figure  4.19  suggests  that  the  mid-range  concentration  of  3x10^5  cm*3  is  nearly 
optimal.  For  higher  concentrations,  the  corresponding  transconductances  at 
all  biases  are  lower,  which  is  due  to  a  combination  of  profile-shape  effects  and 
compensation  effects  on  the  electron  saturated  velocities. 

4.2.3  Diffusion  During  Annealing 

For  Si  and  Se  ions  implanted  to  concentrations  less  than  10^^  cm'^  in 
GaAs.  the  diffusion  constants  for  typical  furnace  anneals  (Ta<900°C)  are  less 
than  10"^'^  cm^/s  (Se)  and  10"^^  cm^/s  (Si)  [60-61].  Ohnuma  et  al.  [62],  using 
capacitance-voltage  techniques  to  measure  electron  profiles  on  Cr-doped 
substrates,  observed  larger  diffusion  constants  for  Si  ions  annealed  under  Si02 
caps.  However,  SIMS  measurements  of  the  Si-atomic  impurity  profiles  in 
undoped  LEC-grown  GaAs  give  much  smaller  diffusion  constants,  suggesting 
that  the  Cr  doping  could  have  influenced  Ohnuma's  results  [63].  Also,  neither 
the  outdiffusion  of  Si  or  Sc  ions  into  Si3N4  or  Si02  caps  nor  the  indiffusion  of  Si 
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from  the  caps  has  been  observed.  Assuming  that  the  diffusion  coefficients 
D(Ta)  for  anneal  temperature  Ta  arc  independent  of  electric-field  and  dopant 
concentration  effects,  the  diffused  concentration  C(x)  after  a  time  t  is 
calculated  from  the  as-implanted  concentration  C'(x)  using  [70] 

C(»)  -  “x-  W)  [^exp[-  .  (4,12) 


This  equation  satisfies  Pick's  diffusion  equation,  with  a  boundary  condition  at 
the  surface  (x=0) 

dx  ■  (4.13) 

corresponding  to  no  transport  of  the  implanted  species  across  the  GaAs/cap 
interface. 

To  first  order,  expressing  C  as  a  Gaussian  with  a  width  a,  the  post- 
anneal  distribution  remains  Gaussian  with  a  width  a,  given  by  [64] 


c7'*  =  a2  +  2Dt.  (4.14) 

This  suggests  that  diffusion  during  annealing  primarily  affects  just  the  width 
of  the  distribution.  Since  MESFET  pinch-off  voltages  arc  roughly  proportional 
to  the  first  moment  of  this  distribution  Rp,  the  larger  width  o  due  to  diffusion 

has  little  effect  on  threshold  voltages,  or  on  other  device  characteristics.  The 
threshold  voltage  begins  to  change  only  when  the  diffusion  length  ~2VDt 
approaches  the  value  of  Rp.  Then,  dopants  diffusing  toward  the  surface  are 
reflected  inward,  whereas  dopants  diffusing  inward  continue  to  move  inward 
and  the  net  moment  of  the  distribution  increases,  leasing  to  larger  pinch-off 
and  more  negative  threshold  voltages. 

Using  this  argument,  we  derive  the  following  criterion  for  how 
diffusion  during  annealing  affects  MESFET  threshold  voltages.  If 

\/4Dt  >  Rp  -  leap  Rps/Rpeap-  (4.15) 

arc  the  projected  ranges  in  amorphous  GaAs  and  in  the 
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Where  Rpj  and  Rpcap 


Where  Rps  and  Rpcap  ^re  the  projected  ranges  in  amorphous  GaAs  and  in  the 
cap,  Tcap  is  the  cap  thickness,  and  Rp  is  the  range  in  crystalline  GaAs  [62].  The 
factor  involving  tcap  describes  shifts  in  the  distributions  due  to  implanting 
through  caps,  as  described  in  ref.  [57].  Figure  4.20  illustrates  this  criterion 


Se  Ions  in  GaAs 


Fig.  4.20.  Threshold  voltages  versus  annealing  time  for  100-,  200- 

and  300-keV  Se  ions  .mplanted  into  bare  wafers  and  for 
300-keV  Se  ions  implanted  through  25-  to  100  nm  Si3N4 
caps.  The  arrows  indicate  the  minimum  time  where  Eq. 
(4.15)  is  satisfied. 

by  showing  calculated  threshold  voltages  versus  anneal  time  for  a  900°C 
anneal  of  implanted  Se.  The  diffusion  constant  used,  10'*'^  cm^/s,  is  for  Se  ions 
at  a  dose  of  approximately  3x10' ^  ions/cm^  [60].  The  arrows  in  Fig.  4.20 
indicate  the  minimum  annealing  time  t  where  Eq.  (4.15)  is  satisfied.  The 
device  threshold  voltages  change  by  less  than  50  mV  up  to  where  Eq.  (4.15)  is 
satisfied,  then  decrease  more  steeply  for  longer  anneals.  Using  this  criterion, 
we  can  conclude  that  diffusion  during  annealing  is  more  important  for  low- 
energy  implants,  and  for  implants  through  thick  caps,  where  the  effective 
moment  of  the  implant  distribution  is  smallest. 
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For  MESFET  technology  the  diffusion  of  Si  and  Se  ions  during  annealing 
is  usually  negligible.  Since  the  diffusion  coefficients  for  Si  and  Se  at  typical 
800-850°C  anneal  temperatures  are  generally  smaller  than  10'***  cm  2/s,  Eq. 
(4.15)  is  rarely  satisfied.  The  use  of  rapid  thermal  annealing  is  becoming 
popular  for  various  reasons.  While  short  anneal  times  do  limit  diffusion,  our 
results  show  that  this  should  not  be  a  major  reason  for  adopting  the  rapid 
thermal  annealing  technique. 

4.2.4  Process  Sensitivity 

4.2. 4.1  Encapsulant  Thickness 
Implanting  through  Si3N4  or  Si02  caps  has  the  beneficial  effect  of 
reducing  the  sensitivity  of  device  threshold  voltages  and  other  characteristics 
to  implant  tilt  and  rotation  angles  [62].  On  the  other  hand,  the  threshold 
voltages  then  become  sensitive  to  the  cap  thickness  and  density.  Figure 
4.21  illustrates  changes  in  threshold  voltage  due  to  changes  in  cap 
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CAP  THICKNESS  x  (nm) 

Fig.  4.21.  Changes  in  threshold  voltages  due  to  changes  in  the  cap 

thickness  for  60-  to  200-keV  Si  ion-implanted  devices  with 
doses  varied  to  achieve  a  constant  -0.5  V  threshold  voltage. 
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thickness  for  MESFETs  implanted  with  60-  to  200-keV  Si  ions  in  the  channel. 
These  calculations  are  for  a  constant  threshold  voltage  of  -0.5  V,  so  for 
increasing  cap  thickness  or  decreasing  implant  energy,  the  ion  dose  is 
increased.  Once  the  dose  is  adjusted,  the  change  in  threshold  voltage  due  to  a 
1-  to  2-nm  change  in  cap  thickness  is  calculated.  Clearly,  the  minimum 
sensitivity  occurs  only  at  vanishing  cap  thickness.  The  sensitivities  are 
largest  for  low  ion  energies,  where  the  ratio  of  cap  thickness  to  the  projected 
range  in  GaAs  is  largest.  To  achieve  a  uniform  wafer-to-wafer  threshold- 
voltage  uniformity  of  <20  mV  for  80-keV  Si  ions  implanted  through  50  nm  caps, 
the  cap  thickness  must  be  controlled  to  within  <lnm. 

4. 2. 4. 2  Recess  Etching 

For  devices  with  recess-etched  gates,  the  etch  depth  must  be  tightly 
controlled  for  good  uniformity.  Figure  4.22  shows  calculated  sensitivities 
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Fig.  4.22.  Changes  in  threshold  voltage  due  to  changes  in  recess  etch 

depth  for  60-  to  200-keV  Si  ion-implanted  devices  with 
doses  varied  to  achieve  a  constant  -0.5  V  threshold  voltage. 
The  implants  are  done  into  unencapsulated  wafers. 
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of  threshold  voltages  due  to  small  etch  depth  variations  for  60-  to  200-keV  Si 
ion-implanted  MESFETs.  Here  the  initial  implants  are  done  into 
unencapsulated  wafers  with  the  doses  adjusted  to  obtain  constant  -0.5  V 
threshold  voltages.  Unlike  the  quadratic  sensitivity  to  cap  thickness,  the 
changes  in  threshold  voltage  with  etch  depth  are  nearly  linear,  and  approach 

constants  when  the  etch  depth  exceeds  the  peak  position  of  the  implant 
concentration  (indicated  by  arrows  in  Fig.  4.22).  Also,  the  sensitivity  to  etch 
depth  is  smaller  than  that  to  cap  thickness. 

4.2.5  Recoil  Ion  Effects 

When  ions  are  implanted  through  Si3N4  and  Si02  layers.  Si,  0  or  N  atoms 
are  knocked  forward  into  the  GaAs.  The  Si  ions  can  be  additional  donors, 
adding  to  the  implanted  donor  concentrations.  The  N  ions  are  electrically 
inactive  in  GaAs,  but  do  contribute  more  damage  in  the  near-surface  region. 
The  O  ions  can  either  form  deep-donor  traps  similar  to  EL2  or  can  complex 

with  and  therefore  compensate  some  of  the  Si  donors,  leading  to  smaller 

electron  concentrations  [65]. 

At  the  shallow  depths  where  the  recoil  atom  concentrations  are  large 
compared  to  those  of  the  implanted  ions,  the  recoil  concentrations  can  be 
predicted  accurately  using  the  TRIM  Monte-Carlo  code  [57,35].  An  extensive  set 
of  TRIM  calculations  were  made  of  recoil-atom  distributions,  resulting  from  Si 
and  Se  implants  through  Si3N4  and  Si02  caps  at  energies  between  60  and  400 
keV  (100  to  400  keV  for  Se).  The  recoil-atom  distributions  were  accurately  fit 
to  the  expression  [66] 

exp(-(x/B)2),  (4.16) 

where  X  is  the  oepiii  ui  ♦he  G2.\s  in  nm,  and  n  (--  1  to  1.5),  B  and  A  are 
parameters  dependent  on  the  ion  energy,  the  cap  thickness,  the  implant  ion, 
and  the  recoil  atom  (Si,0,  or  N). 

Figure  4.23  shows  calculations  of  Si  recoil-atom  concentrations 

resulting  from  a  300-keV  Se  implant  through  a  100  nm  Si3N4  overlayer.  The 
parameter  B  for  this  implant  is  approximately  300  nm,  therefore  the  recoil- 
atom  concentrations  fall  off  as  a  power  law  function  of  depth  for  most  depths 
shown  in  Fig.  4.23. 

Figures  4.24  and  4.25  show  the  effects  of  Si  recoil  atoms  on  MESFET 
threshold  voltages.  For  these  calculations,  the  ion  dose  was  adjusted  to  obtain 
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an  initial  threshold  voltage  of  -0.5  V,  then  the  recoil  atom  distribution  was 
added,  and  the  threshold  voltage  was  recomputed.  For  a  Si  ion  primary 
implant,  the  threshold  voltages  become  as  much  as  60  mV  more  negative,  while 
for  the  larger-mass  Se  ions,  300  mV  changes  arc  found.  For  a  Be  implant,  with 


DEPTH  (pm) 

Fig.  4.23.  Concentrations  of  Se  ions  implanted  through  100  nm  Si3N4 

caps  and  the  concentrations  of  knocked-on  Si  atoms  from 
the  cap. 

a  typical  buried-p-layer  implant  dose,  the  change  in  threshold  voltage  due  to 
recoils  is  less  than  3  mV.  The  maximum  cap  thicknesses  tcap  considered  in 
these  calculations  are  about  l.SRps,  where  Rps  is  the  projected  ion  range  in 
amorphous  GaAs  [57]. 

4.2.6  Conclusions 

In  this  paper  we  have  addressed  several  ion-implantation  effects  on 
GaAs  MESFET's.  With  decreasing  channel  implant  energy  E,  transconductances 
increase  less  quickly  than  1/E,  the  increase  over  a  factor  of  10  reduction  in  E  is 
found  to  be  only  a  factor  of  3  to  4.  At  low  gate  biases,  transconductances  for 
implanted  devices  ar>.  smaller  than  for  box  profiles.  The  presence  of  substrate 
acceptors  or  implanted  buried  p  layers  is  found  to  improve  low-bias 
transconductances. 
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A  criterion  has  been  derived  to  determine  whether  or  not  dopant 
diffusion  during  post-implant  annealing  affects  device  characteristics.  For  Si 
and  Se  at  typical  channel  implant  doses  and  energies,  and  with  typical  800- 
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Fig.  4.24.  The  effects  of  knocked  -on  Si  atoms  on  device  threshold 

voltages  when  60-  to  200-keV  Si  ions  are  implanted 
through  Si3N4  caps.  Without  recoils,  a  -0.5  V  threshold 
voltage  was  obtained  for  all  ion  energies  and  cap 
thicknesses,  then  the  recoil  distributions  were  added, 
resulting  in  the  threshold  voltages  shown. 


0  50  100  150 

CAP  THICKNESS  (nm) 


Fig.  4.25. 


Same  as  Fig.  4.24  for  Se  ion  implantation  through  Si3N4 
caps. 


850°C  furnace  or  rapid  thermal  anneals,  dopant  diffusion  effects  on  MESFET 
threshold  voltages  are  insignificant. 

Calculations  are  presented  for  the  effects  of  cap-thickness  and  recess- 
etch  depth  fluctuations  on  threshold  voltages.  Also,  for  implantation  through 
these  Si02  and  Si3N4  layers.  Si  atoms  from  the  cap  which  arc  knocked  forward 
into  the  GaAs  change  threshold  voltages  by  up  to  60  mV  for  Si  or  300  mV  for  Se 
ion  implantation  (for  a  target  threshold  voltage  of  -0.5  V). 

A  quantitative  understanding  of  these  effects  should  greatly  benefit  the 
design  and  optimization  of  GaAs  MESFETs  for  digital  and  monolithic-microwave 
applications. 
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V. 


DEVELOPMENT  OF  PROCESS  MODEL 


5.1  Background 

There  is  currently  a  need  for  process  and  device  modeling  tools  in  the 
GaAs  monolithic-microwave  integrated-circuits  industry.  Crucial  to  the 
evolution  of  high  yield  GaAs  MMIC’s  the  development  of  processes  that 
tightly  control  device  electrical  characteristics,  such  as  ion-implanted 
MESFET  threshold  voltages  or  saturated  currents.  Process  and  device 
modeling  can  provide  a  quantitative  understanding  of  materials  effects 
(wafer  impurities,  dslocations),  ion-implant  effects  (wafer  tilt  and  rotation 
angles,  annealing  encapsulant  effects),  or  how  day-to-day  process  variations 
affect  circuit  uniformity.  Also,  the  design  of  MMIC's  is  crucially  dependent 
on  device  modeling  tools  that  can  accurately  predict  circuit  performance 
before  fabrication.  Traditionally  MESFET's  are  fabricated  then  measured  and 
the  device  characteristics  are  fltted  to  mathematical  expressions,  suitable  for 
use  in  SPICE  [67-69]  or  Harmonic  balance  [70]  codes  to  predict  circuit 

performance.  As  the  measurement  techniques  used  to  obtain  small-signal 

parameters  are  not  always  unambiguous,  device  modeling  codes  that  predict 

these  quantities  are  of  great  benefit. 

In  this  section  we  describe  a  combined  process  and  device  model 

developed  under  this  program  for  GaAs  MESFET  technology:  Instead  of 

attempting  to  solve  ab  initio  equations  for  the  transport  of  ions  in  matter  for 
ion  implant  process  modeling  or  multi-dimensional  Poisson's  and  current- 
continuity  equations  for  device  modeling,  our  code  uses  analytical  models. 
More  than  10  years  of  GaAs  materials  and  device  research  has  resulted  in  the 
availability  of  many  useful  process  and  device  models  for  GaAs  MESFETs.  The 
emphasis  in  the  development  of  the  model  is  the  integration  of  these  into  a 
software  package  that  provides  an  easy-to-use  tool  for  non-modeling- 
professional  process  development  or  manufacturing  engineers. 

Figure  5.1  we  show  an  overview  of  the  model  [71],  The  process  model 
section  calculates  ion-implant  impurity  profiles  in  one  dimension,  and 
includes  the  diffusion  of  dopants  during  annealing,  recess  etching,  and  the 
indiffusion  of  dopants  from  the  surface  (used  in  some  JFET  fabrication 
processes).  To  model  materials  effects,  background  shallow  acceptor  and 
donor  [72],  EL2  [60],  and  deep-acceptor  (e.g.  Cr  [20])  concentration  profiles 
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Fig.  5.1.  Overview  of  GATES. 

can  be  constructed.  Piezoelectric  effects  are  calculated  using  the  dielectric- 
edge  approximation  [73],  and  short-channel  effects,  due  to  the  encroachment 
of  n'*'  implant  profiles  into  the  channel  region,  are  calculated  using  a  method 
proposed  by  Chen  et  al.  [74].  Threshold  voltages  and  carrier  profiles  are 
obtained  from  either  a  Poisson  solver  [8]  or  an  abrupt-depletion 
approximation  for  arbitrary  carrier  profiles  [5],  and  Ids  is  calculated  using 
the  two-piece  mobility  approximation  [58,75],  where  compensation  effects  on 
mobility  and  saturated  velocity  are  included  [75,76].  For  the  source  and  drain 
resistances  for  four  general  types  of  MESFET  topologies,  we  calculate  layer 
resistivities  for  arbitrary  implants,  then  we  sum  the  resistances  of  the 
various  regions  between  the  contact  and  the  gate,  according  to  the  specific 
FET  geometry.  The  contact  resistance  is  calculated  using  a  ND"^  dependence 
[77],  though  usually  the  user  must  overide  the  default  values,  due  to  the  large 
process  variations  in  Au/Ge  contact  resistances  [77].  Also,  all  of  the 
transconductanccs,  capacitances,  and  resistances,  used  to  determine  small-  or 
large-signal  characteristics  arc  obtained  analytically.  For  parameter 
extraction,  the  I-V  curves  can  be  fit  to  SPICE  formulas  proposed  by  Curtice 
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[67],  Staiz  et  al.[68],  or  Larson  [69].  Concentration  profiles  or  electrical 
characteristics  can  be  plotted. 

Sections  5.2  and  5.3  describe  the  process  and  device  models 
incorporated  into  the  code.  Piezoelectric  and  short-channel  effects  are  not 
discussed,  as  the  methods  used  are  similiar  to  those  of  Chen  et  al.  [74].  As  with 
any  modeling  program,  calibrating  some  physical  parameters  is  of  vital 
importance. 

5.2  Process  Models 

5.2.1  Ion  Implantation 

The  code  allows  for  four  different  kinds  of  profiles  in  one  dimension; 
implants  for  the  channel  donor  Nd‘,  the  source  and  drain  regions  Nd"*’, 
another  donor  implant  No',  and  an  acceptor  implant  Na-  Multiple  implants 
can  be  added  to  any  one  of  these  profiles. 

Pearson-lV  functions  [47]  are  used  to  model  the  profiles  of  Be,  Si,  and 
Sc  ions,  [See  Section  4.1].  For  maximum  cross  wafer  device  uniformity  and 
the  narrowest  profiles  for  a  given  ion  energy,  the  code  assumes  implants  are 
performed  at  a  tilt  angle  of  at  least  9°  and  a  rotation  angle  of  45°  from  the 
major  flat  (Oil). 

In  calculating  Be  profiles  in  GaAs-implanted  through  Si02  or  Si3N4  caps, 
measured  profiles  for  implantation  into  bare  wafers  are  shifted  toward  the 
surface  by  an  amount  tRps/RpC.  where  t  is  the  cap  thickness,  and  RpS  and 
RpC  are  the  projected  ranges  in  the  semiconductor  and  cap  [57]  (Fig.  5.2). 


Fig.  5.2, 


Measured  concentration  profiles  for  50-keV  Be  ions  implanted 
into  bare  GaAs  wafers  and  into  wafers  with  45  and  79  nm  CVD 
Si3N4  caps,  compared  with  fits  using  the  range-scaling  theory. 


Moments  of  ion  distributions  in  the  amorphous  cap  are  interpolated  from 
tables  of  values,  calculated  using  the  TRIM  Monte-Carlo  code  [35].  The  dose  in 
the  GaAs  is  normalized  to  incident  ion  dose  reduced  by  that  lost  in  the  cap. 
Consistent  with  TRIM  results,  a  Pearson-I  distribution  [47]  is  asssumed  in  the 
cap. 

For  Si  and  Se  ions,  implanting  through  caps  reduces  the  magnitudes  of 
the  channeling  tails.  This  effect  is  simulated  by  mixing  a  shifted  Pearson-IV 
function  with  a  shifted  TRlM-calculated  Gaussian  function,  using 

C(x)  =  (I-p)CGaus(y)  +  PCpiv<y).  (5  1) 

where  y=x-tRps/RpC’  Cpiv  are  normalized  to  1  at  x=Rp.  C(x)  is 

normalized  to  the  dose  in  the  GaAs.  This  method  is  illustrated  in  Fig.  5.3.  The 
mixing  factor  p  (<1.)  depends  on  the  ratio  of  the  cap  thickness  t  to  Rpg.  but 
has  only  been  characterized  for  a  few  cases  [57]. 


Fig.  5.3.  Measured  concentration  profiles  for  100-keV  Si  ions 

implanted  into  GaAs  through  a  79nm  CVD  Si3N4  cap, 
compared  with  fits  using  Eq.  (5.1). 

When  ions  are  implanted  through  Si02  and  Si3N4  caps.  Si,  0,  or  N  ions 
are  knocked  into  the  GaAs,  where  the  Si  ions  are  additional  donors, 
increasing  carriei  densities.  Knock-on  atom  distributions  arc  calculated 
iisi.ig  TRIM  [35],  and  arc  fit  to  the  expression 

Cj(x)  =  —  exp(-(x/a)*), 

xD  +  b  (5.2) 
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where  b=0.1  and  x  is  the  depth  in  nm.  The  parameters  A, a,  and  n  for  Si,  0,  or 
N  atoms  can  be  interpolated  from  tables  as  a  function  of  cap  (Si02  or  Si3N4), 
cap  thickness  (0<t<1.3RpS),  ion  (Si  or  Se),  and  ion  energy  (60  to  400  keV).  The 
knock-on  concentrations  arc  only  significant  compared  to  the  implant-ion 
distribution  near  the  surface  (where  channeling  tails  arc  unimportant),  and 
MESFET  pinchoff  voltages  are  increased  by  <1%  (Be),  5-10%  (Si),  and  20-30% 
(Se). 

5.2.2  Substrate  Materials 

Four  significant  impurites  can  be  found  in  liquid-encapsulated 
Czochralski  (LEC)  or  horizontal-Bridgman  (HB)  grown  GaAs  substrates 
[72,74,78],  Table  5.1  gives  typical  ranges  and  default  values  of  these 


Table  5.1.  Wafer  Impurity  Concentrations  (10^^  cm"3). 


LEC  wafers 

HB  wafers 

Impurity 

Ion 

Range 

Default  Range 

Default 

Shallow  Acceptor 

C 

1-6 

3 

«Si 

0 

Shallow  Donor 

Si 

«C 

0 

5-15 

10 

Deep  Donor 

EL2 

5-20 

20 

«Cr 

0 

Deep  Acceptor 

Cr 

«EL2 

0 

20-100  50 

in  modem,  manufactured 

wafers. 

Both 

types  of 

wafers  are  : 

In  LEC  wafers. 

the 

Fermi 

level  is 

pinned 

mid-gap  by 

insulating. 

compensation  of  the  dominant  shallow  acceptor  by  the  native  deep-donor 
defect,  EL2.  In  HB  wafers,  the  dominant  donor  impurities  arc  compensated  by 
added  mid-gap  acceptor,  Cr.  The  code  allows  varying  these  concentrations,  to 
observe  their  effect. 

5.2.3  Diffusion 

Dopant  diffusion  during  post-implant  annealing  is  calculated  using 

[47] 


C(x)  = 


-oo 

- 

dx'  C(x') 

Jo 

exp 

- 

+exp 

(x+x')*! 

4Dt 

4Dt  J 

(5.3) 


where  C  and  C  are  the  as-implanted  and  annealed  distributions,  D(T)  is  the 
assumed  constant,  concentration-independent  dffusion  coefficient  at  the 
anneal  temperature  T,  and  t  is  the  anneal  time.  This  expression  is  a  solution 


of  Pick's  equation  when  no  outdiffusion  from  the  substrate  is  present 
(dC/dx=0  at  x=0). 

The  diffusion  constants  for  low-concentration  Si,  Se  and  Be  ions  in 
GaAs  are  so  small  that  dopant  diffusion  during  typical  furnace  anneals  is  an 
insignificant  factor  in  MESFET  fabrication.  MESFET  threshold  voltages  are 
affected  when  the  dffusion  length  V4Dt  exceeds  the  mean  ion  range  in  the 
semiconductor,  RpS'^RpS/RpC  [79].  With  typical  diffusion  constants  for  Se 
ions  at  channel  concentrations  of  10*7  cm‘3,  (10'*^  cm^/s  at  900“C  [60]), 
anneal  times  greater  than  20  min  for  bare-wafer  implants  are  needed  to 
affect  threshold  voltages  more  than  a  few  percent.  Si  diffusion  constants  are 
not  well  characterized,  but  are  much  smaller  [61,63,80].  Buried  p  layer 
diffusion  could  also  affect  threshold  voltages,  but  diffusion  constants  for  Be 
or  Mg  ions  at  concentrations  <10*^  cm'^  are  not  well  characterized  [81,23,82]. 
In  one  case,  SIMS  measurements  detected  no  diffusion  in  a  Be  buried  p  layer 
for  a  typical  furnace  anneal  [23].  It  is  not  well  known  whether  the  constant- 
D  assumption,  implicit  in  Eq.(5.3),  is  valid  at  these  low  dopant  concentrations, 
but  this  assumption  has  generally  been  used  in  the  extraction  of  reported 
diffusion  coefficients. 

5.2.4  Recess  Etching 

Recess  etching  is  modeled  by  shifting  the  profiles  toward  the  surface. 
Automatic  etching  to  a  desired  threshold  voltage  or  saturated  current  can  be 
done. 

5.2.5  Epitaxial  Layers 

Competing  with  ion-implanted  MESFET  technology  is  epitaxial-doped- 
layer  growth,  best  controlled  by  molecular  beam  epitaxy  or  MOCVD  methods. 
Layer  thicknesses  and  doping  concentrations  can  be  input  to  the  code. 
Implant  profiles  for  the  source  and  drain  regions  can  then  be  added  to  these 
profiles  to  simulate  epitaxial-fabricated  MESFET's. 

5.2.6  Activation 

Typical  Si  and  Se  donor  activation  coefficients  are  between  70%  and 
957o,  and  at  the  small  concentrations  used  for  buried  p  layers,  near  100% 
acceptor  activation  is  usually  found.  The  default  donor  and  acceptor 
activations  used  arc  therefore  80%  and  100%  for  Si  and  Se,  respectively.. 
Within  the  expected  range  of  possible  activation  achieved  with  different 
annealing  caps  [82],  temperature  and  dose  [83],  and  other  factors,  greatly 
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varying  threshold  voltages  are  predicted.  The  donor  activation  can  be 
calibrated  using  measured  threshold  voltages. 

For  the  high  donor  concentrations,  used  for  Nd"*"  implants  and 
sometimes  for  low-energy,  high-dose  channel  implants,  the  activation  may 
be  concentration  or  depth  dependent  [84].  At  low  concentrations,  the 
electron  and  donor  concentrations  differ  by  a  constant  factor,  but  above  a 
certain  donor  concentration,  the  electron  profile  is  nearly  flat  (Fig.  5.4.). 


DEPTH  (gm)' 

Fig.  5.4,  Implant  profile  for  200-keV  Si  ions  implanted  with  a  dose 

of  4x10*'^  ions/cm^,  and  GATES  calculations  of  electron 
carrier  densities  for  anneal  temperatures  of  800°,  900°, 
and  1000°C,  using  Eq.  (5.4).  The  reason  for  the 
difference  at  large  depths  is  not  known. 


This  saturation  is  modeled  by  assuming  that  for  amphoteric  Si  ions,  high- 
concentration  activation  is  controlled  by  the  distribution  of  Si  between  Ga 
and  As  sites  Assuming  equilibrium  conditions,  this  is  calculated  using  [85,13] 


(S'As  1 


_1_ 

K. 


SL 

4 

’h2 

"i 

(5.4) 


[SiGa"^! 

where  Ky  is  the  equilibrium  constant  for  the  reaction  between  Si  on  As  and 
Ga  sites  and  Ga  and  As  vacancies  (sec  calibration  in  Sect.  IV),  [Si]  is  the  total 
Si  concentration,  and  nj  is  the  intrinsic  electron  concentration  at  the  anneal 
temperature.  This  model  predicts  that  y  increases  from  0  as  [Si]  increases. 


therefore  the  activation  coefficient  "H  =(  1 -Y)/(  1 +Y)  decreases  from  unity. 

Since  the  low-concentration  activation  coefficient  TJo  is  less  than  unity,  we 
hypothesize  a  minimum  Y=(l'Tlo/(l+'no)-  The  model  predicts  that  the  saturated 
electron  concentration  increases  (therefore  n"*"  layer  resistivities  decrease 
[86])  with  higher  anneal  temperature.  (The  opposite  effect  can  occur  is  As 

loss  is  significant.  Then  Ky  decreases  with  higher  temperature,  leading  to 

larger  Y  values,  and  lower  activation.) 

5.3  Device  Modeling 

5.3.1  Threshold  Voltages 

To  simulate  device  I-V  characteristics,  carrier  densities  are  required, 
which  are  calculsted  from  the  derived  impurity  profiles  by  solving  Poisson's 
equation  in  one  dimension  using  a  Newton  iteration  method  [8],  and  using 
the  boundary  condition  at  the  surface,  dn(x)/dx=0.  This  gives  a  gate-bias- 
independent  electron  profile  resembling  the  active  donor  profile  (Fig.  5.), 
except  near  the  surface,  where  the  boundary  condition  forces  n(x)  to  depart 
from  the  usually  steeply  increasing  profile,  and  at  large  depths,  where  the 

substrate  junction  acts  to  cut  off  the  tails  in  the  implant  profiles  (Fig.  5).  The 
MESFET  pinchoff  Vp  and  threshold  voltages  Vth  are  defined  as  [5] 


xdx  „(x),  Vai-Vy-Vp, 


(5.5) 


where  Viji-O.TS  V  is  the  Schottky  barrier  height,  which  is  nearly  constant 
over  widely  varying  process  and  operating  conditions  [86] 

5.3.2  Layer  Resistivity 

The  conductivity  of  an  implanted  layer  is  given  by  [86] 


<7  =  qW/L 


I 


oo 


dx  n(x)  /i(x). 


(5.6) 


where  W  and  L  are  the  layer  width  and  length,  and  x)  is  the  electron 

mobility,  which  depends  on  the  carrier  concentration  and  on  the 
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Fig.  5.5.  Comparison  of  measured  donor  profiles  (solid  line), 

resulting  calculated  electron  profiles  (solutions  to 
Poisson’s  equation  for  n(x)  with  zero  bias;  dashed  line), 
and  measured  electron  profiles,  determined  using  the  CV 
method  (points).  Also  shown  are  the  EL2,  C,  and  Cr 
concentrations  used  in  the  calculations  of  n(x). 


compensation  ratio  at  each  depth  [75,76].  Similiarly,  we  define  the  saturated 
current  as 


f"’ 

^sat  =  dW  dx  n(x)v5(x), 

Jx„  (5.7) 


where  Vg  is  the  saturated  electron  velocity,  also  depending  on  the  carrier 
density  and  compensation,  therefore  the  depth.  The  code  uses  Caughey- 
Thomas  expressions  for  electron  mobilities  [87]  which  give  good  fits  to  the 
values  tabulated  by  Walukiewicz  et  al.  [76].  The  electron  saturated  velocity  is 
related  to  the  mobility  using  Vsat=(0.64+ 1 .3  8p. -0. 82n 2)x  i  o7cm/s  ([1  is  in 

m  2/v/s),  which  was  suggested  by  Xu  and  Shur  [88].  The  constants  were 
adjusted  to  obtain  slightly  better  agreement  with  Xu  and  Shur's  calculations 
over  the  small  range  of  mobililties  normally  encountered. 

6  1 


GaAs  surfaces  are  pinned  by  the  presence  of  a  surface  potential  Vgf, 
which  is  approximately  0.6  V,  independent  of  most  types  of  layers  on  the 
surface  [89].  This  surface  potential  depletes  carriers  immediately  below  the 
surface,  so  that  instead  of  integrating  from  x=0  in  Eqs.(5.6)-(5.7),  the  lower 
limit  is  Xq,  which  in  the  abrupt  depletion  approximation  [S],  is  obtained  from 


r*o 

Vjf  =  xdxn(x), 

^o^r  Jo 


(5.8) 


5.4  Conclusions 

We  have  described  an  approach  to  modeling  GaAs  MESFET's  that  makes 
extensive  use  of  analytical  models.  This  approach  provides  a  tool  which  has 
allowed  calculations  of  parameters  useful  in  correlation  substrate  and 
process  defects  with  device  characteristics. 
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V I .  COMPARISON  OF  THRESHOLD  VOLTAGE  FLUCTUATIONS  FOR  Si 

AND  Se  IMPLANTED  GaAs  MESFETs 

6.1  Background 

The  effect  of  process  variables  on  LSI  circuits  is  an  unresolved  issue  in 
GaAs  technology.  The  final  parameters  arc  affected  by  each  of  the  process  and 
material  variables.  The  standard  deviation  of  any  MESFET  parameter  across  a  3 
inch  wafer  will  vary  according  to  the  summation  of  the  standard  deviation  of 
each  variable.  LSI  design  requires  a  limited  variance  on  the  parameters  of 
these  FETs  in  order  to  achieve  functionality  of  the  desired  circuits.  The 
outcome  of  this  scenario  is  that  a  manufacturing  line  establishes  a  process  and 
limits,  as  much  as  possible,  any  changes  to  that  process. 

In  the  design  of  an  experiment  one  must  choose  which  steps  of  the 
process  can  be  va  ied  without  disturbing  the  process  or  equipment.  Thus,  we 
chose  the  active  implant  species  as  one  process  variable  in  our  experiment. 
Both  silicon  and  selenium  were  chosen  as  channel  implant  species,  since  both 
are  routinely  used  at  the  Rockwell  GaAs  pilot  line  and  should  give  different 
results  in  the  planned  experiments. 

One  topic  disputed  in  the  literature  is  the  effect  of  dislocation  density  on 
the  threshold  voltages  of  MESFETs  fabricated  on  GaAs.  Another  similar  issue  is 
EL2  (a  deep  level  trap)  which  correlates  after  crystal  growth  with  dislocation 

density.  That  is,  dislocation  density  and  EL2  both  demonstrate  a  "W"  pattern  in 
distribution  on  LEG  grown  wafers  throughout  a  boule.  The  EL2  distribution, 
however,  will  level  off  following  a  long  term  anneal  of  10  hours  at  850C.  It  has 
been  assumed  that  the  anneal  after  implant  will  also  diffuse  the  distribution. 

Industry  wide  whole  boule  annealing  was  initially  initiated  to  reduce  residual 
stress  in  boulcs  to  limit  cracking  during  grinding  and  slicing.  The  Rockwell 

grown  material  used  in  these  experiments  was  grown  under  reduced 
temperature  gradients,  and  thereby  there  was  no  concern  of  cracking  and  it 

therefore,  did  not  receive  a  WHOle  Boule  Anneal  (WHOBA).  No  boule  grown 

under  a  heat  shield  at  Rockwell  has  ever  cracked  during  processing.  Thus, 
that  these  wafers  had  not  been  whole  boule  annealed  added  an  unforeseen 
material  variable  not  realized  until  the  final  stages  of  the  effort. 

Another  material  variable  over  which  there  is  little  control  is 

stoichiometry,  the  atomic  ratio  between  Ga  and  As  in  the  boule.  This  is  a 

known  contributor  to  EL2,  however,  the  control  of  stoichiometry  has  always 
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been  done  on  a  macroscopic  level  with  LEC  crystal  growth.  The  growth  runs 
are  monitored  for  loss  of  arsenic  via  loss  in  gross  weight  of  the  charge  which 
is  compensated  for  by  extra  arsenic  being  added  to  the  initial  charge.  As  can 
be  seen  from  the  As-Ga  phase  diagram  this  will  not  control  stoichiometry  on  a 
macroscopic  level.  What  is  needed  is  to  control  the  As  pressure  directly  with  is 
a  separate  As  source.  Industry's  approach  to  the  dislocation  question  has  been 
to  dope  GaAs  with  In  to  lower  the  dislocation  density.  This  works  by  the  larger 
In  atom  substituting  for  Ga  and  blocking  the  propagation  of  dislocations 
through  the  material. 

Thus  our  experiments  contained  one  process  variable,  implant  species, 
and  three  material  variables:  In-doping,  whole  boule  annealing,  and  growth 
technique  (LEC  vs  VGF).  The  first  lot  contained  4  wafers  of  non-whole 
annealed  In-doped  low  dislocation  material  and  4  wafers  of  annealed  high 
dislocation  standard  material.  The  implant  species  for  this  lot  was  Se.  Lot  #2 
was  identical  to  lot  #1  except  the  active  channel  was  implanted  with  Si.  Lot  #3 
would  contain; 

Lot  it3  (unfinished) 

4-  standard  high  dislocation,  WHOBA  LEC  wafers; 

2-  In  doped,  low  dislocation,  non-WHOBA  LEC  wafers; 

2-  In  doped,  low  dislocation,  WHOBA  LEC  wafers; 

2  -  non-doped,  low  dislocation,  WHOBA  VGF  wafers. 

This  lot  was  not  completed  due  to  exhaustion  of  the  processing  and 
measurement  funds  by  lots  1  and  2. 

6.2  Physical  Modeling 
6.2.1  Model 

In  the  interim  technical  report  issued  for  this  contract  (90),  we 
proposed  a  quantitative  relation  for  carrier  activation  in  Si  ion  implanted  GaAs 
based  on  results  obtained  from  several  authors  (91).  This  model  assumed  that  a 
localized  defect  related  to  stoichiometry  imbalance,  namely  Ga  vacancies 
('•^Ga)-  was  responsible  for  causing  variations  observed  in  the  activation  of  Si 
implanted  GaAs.  The  value  for  the  activation  ratio  as  a  function  of  Vq  a 
obtained  was 


1  -  A/Vca 
'I  +  A/VCa 


(6.1) 
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where  A  is  a  material  and  temperature  related  constant.  Using  the  data  of  (91) 
relating  MESFET  threshold  voltage  to  the  proximity  from  a  dislocation  for  Si 
channel  implants,  we  derived  a  VQa  concentration  dependence  as  a  function  of 
the  distance  r  to  a  dislocation, 

VCa/A  =  Ab  +Ao  exp  [-(r/a)J-5]  (6.2) 

In  Eqn.  6.2,  Ab  is  the  background  ratio  of  Ga  vacancies  in  the  material  and  Aq 
and  a  are  chosen  for  a  best  fit  of  the  experimental  data.  The  background  ratio, 
Ab,  is  chosen  to  give  the  correct  average  percent  carrier  activation  and 

depends  upon  the  process  and  implanted  atom.  For  instance,  to  achieve  a  70% 
activation  for  Si  implanted  GaAs  annealed  under  a  Si3N4  cap,  Ab  =  5.6. 

The  activation  of  Si,  a  Group  IV  dopant,  in  GaAs  is  strongly  dependent  on 
the  availability  of  Vca-  For  wafers  in  which  there  is  a  non  uniform  variation 
in  VGa,  a  resultant  variation  in  the  activation  of  electrically  active  Si  should 
result.  Selenium,  however,  being  a  Group  VI  element  should  be  much  less 
dependent  on  the  local  concentration  of  V^a,  resulting  in  a  small  variation  in 
activation  when  compared  to  that  of  the  silicon.  The  net  result  should  be  an 

increased  fluctuation  in  MESFET  electrical  parameters  for  wafers  in  which  a 

nonuniform  distribution  of  VQa  exists.  In  the  following  we  describe  a 
numerical  and  experimental  approach  with  which  we  verify  this  hypothesis. 

6.2.2  Numerical  Calculations 

Before  the  Monte-Carlo  calculations  can  be  made  we  must  estimate  the 
base  activation  of  the  implanted  channel  layers.  We  use  our  process  model 

GATES  to  derive  the  base  activation  of  the  Si  and  Se  ions  obtained  from  actual 
measurements  on  the  devices.  For  Se,  an  implantation  of  320  keV  ions  at  a  total 
dose  of  2.2x10^2  cm'2  through  a  750A  Si3N4  cap  of  density  3  3  is  used. 

Including  the  Si  recoils  from  the  cap  in  the  total  active  dose,  an  80%  activation 
is  required  to  achieve  the  experimentally  observed  average  threshold  voltage 
of  Vj  =  -0.55  V.  Although  piezoelectric  effects  are  not  completely  negligible, 
even  for  Lg  =  2  mm,  we  neglect  these  in  the  present  calculations.  For  the  Si 
implants,  using  145  keV,  identical  cncapsulant  and  a  dose  of  2.2x10^2  cm’2.  we 
arrive  at  an  activation  of  64%  to  fit  the  experimentally  observed  Vt  of  -0.4  V. 

Equation  6.2  forms  the  basis  for  our  Monte-Carlo  calculations.  Although 
we  do  not  expect  the  parameters  derived  from  the  experimental  data  u.sed  to  be 
universally  valid,  they  are  used  for  all  the  calculations  performed  following. 
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Also,  since  they  are  likely  to  vary  in  unknown  ways  with  substrate  factors 
(e.g.  Ga-As  stoichiometry,  background  impurities,  boule  annealing,  etc.)  and 
processing  (type  of  anneal,  encapsulant,  anneal  time,  ion  dose  and  energy, 
etc.)  we  expect  only  qualitative  agreement  with  measurements. 

We  model  the  device  to  dislocation  distance  as  a  random  ’.ariable.  This  is 
justified  since  dislocations  are,  to  first  order,  randomly  distributed.  We  assume 

that  the  dislocations  are  distributed  randomly  over  a  wafer  area  of  dimensions 
y  =  200  pm  by  Np,  where  Np  is  the  number  of  FET's  averaged.  The  FETs  are 
arranged  in  a  row  at  y=100pm  and  are  a  distance  p  apart.  For  any  FET  in  the 
row,  Voa/A  is  obtained  by  summing  Eqn.  6.2  over  the  randomly  distributed 

dislocations  in  the  area  which  have  r  <  100  |Im.  One  then  obtains  the  spatially 

dependent  activation  T)  for  the  FET  and  calculates  its  threshold  voltage.  A 

simple  average  of  the  threshold  voltages  of  the  Np  FET's  is  then  calculated  and 
the  standard  deviation,  AVp  obtained.  Figure  6.1  shows  the  threshold 


Fig.  6,1.  Plot  of  threshold  voltage  fluctuation,  AVj,  versus  EPD  for 

Si  and  Se  implanted  GaAs  MESFETs.  The  plots  arc  a  result  of 
a  Monte-Carlo  calculation  which  takes  into  account  the 
effect  of  dislocations  on  MESFET  devices.  The  background 
activations  assumed  for  the  calculations  were  64  and  80% 
for  Si  and  Se,  respectively. 


voltage  standard  deviation  for  500  FETs,  placed  p  =  10  pm  apart  for  Si  and  Sc 
channel  implants  versus  Etch  Pit  Density  (EPD).  For  the  Si  implanted  wafers 
with  T|o  =  64%,  Vt  decreases  with  increasing  EPD,  and  AVp  varies  from  zero  at 
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zero  EPD  (no  dislocations  gives  no  threshold  voltage  fluctuations  in  the 
present  model),  to  a  maximum  at  EPD  =  10^  to  2x10^  cm'^,  falling  off  rapidly  at 
higher  EPD.  For  the  Se  implanted  wafers,  a  value  of  T^o  =  80%  is  used  and  it  is 

seen  that  the  AV-f  variations  are  small  (<0.01  V)  at  EPDs  =  10^  cm‘2  gradually 
increasing  to  =  0.015  V  at  EPDs  =  5x10^  cm"^. 

6.3  Experimental  Approach 

6.3.1  Test  Mask  Design 

The  mask  was  laid  out  with  a  line  of  MESFETs  as  dense  as  the  Rockwell 
depletion  mode  design  rules  would  allow.  A  group  of  ten  FETs  arranged  with 
common  gate  pads,  common  source  pads  and  common  backgates  with  separate 
drains  were  grouped  together  in  a  dense  row  pattern  (DRP).  The  MESFETs  were 
5  p  wide,  3  p  in  gate  length  separated  by  10  p.  Three  rows  of  80  DRPs  were  laid 
in  a  field.  There  are  44  fields  per  wafer.  Therefor  a  single  row  across  a  3  inch 
wafer  crosses  8  fields  or  640  DRPs  or  6400  FETs.  Five  different  backgate  lengths 
were  split  among  the  fields  with  backgate  lengths  of  5,  10,  15,  20,  and  50  p 
6.32  Process  Description 

The  process  chosen  for  the  fabrication  of  the  MESFETs  was  the  then 

standard  depletion  mode  process  at  the  Rockwell  pilot  line.  Prior  to  processing, 
incoming  material  is  qualified  and  characterized.  This  involves  using  a 

Schottky  diode  mask  set  on  two  representative  wafers  from  each  boule.  This 

yields  information  about  the  activation  of  the  implant  in  this  material.  The 

boulcs  arc  rejected  or  accepted  based  on  C-V  measurement  results  obtained 
from  these  diodes. 

A)  After  wafer  qualification  and  acceptance,  the  GaAs  surface  is 
prepared  and  a  Si3N4  cap  is  deposited.  The  implant  profile  is  dependent  on  the 
Si3N4  cap  thickness  and  therefore  MESFET  threshold  uniformity  can  be 
affected  by  the  reproducibility  of  the  Si3N4  cap  thickness.  Fig  6.2  shows  the 
excellent  control  over  Si3N4  thickness  at  the  pilot  line.  This  figure  shows  the 
variation  in  cap  thickness  over  41  wafer  lots.  Each  data  point  represents  a  lot 

of  5  to  20  wafers.  After  cap  deposition,  implant  regions  are  defined  and 
implanted.  Two  N-implants  arc  required.  The  first  being  the  standard 

selenium  implant  for  the  channel  of  the  FETs  and,  as  stated  earlier,  was  the 
process  variable  of  the  experiment,  the  exchange  of  silicon  for  selenium.  The 

second  N-  implant  is  followed  by  a  deep  p-well  implant  through  the  N-  implant 

mask.  The  p-wcll  implant  is  necessary  to  improve  the  radiation  hardness  of  the 
array . 
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Fig.  6.2.  Sputtered  Si3N4  cap  thickness  for  41  lots. 

B)  After  implant  definition,  an  Si02  layer  is  deposited  using 
plasma  enhanced  chemical  vapor  deposition  (PECVD).  The  Si02  layer  prevents 
degradation  of  the  GaAs  during  implant  anneal.  After  Si02  deposition,  the 
wafers  are  annealed  to  remove  the  implant  damage. 

C)  Ohmic  and  Schottky  contacts  are  made  after  anneal.  The 

contacts  are  defined  in  resist,  the  Si02  cap  is  opened  up  using  reactive  ion 
etching  (RIE)  and  the  Si3N4  cap  is  opened  up  using  a  barrel  plasma  process. 
After  opening  the  dielectric,  ohmic  metal  is  deposited  and  lifted  off.  The 
dielectrics  assist  in  the  lift-off  process  by  providing  an  undercut  profile. 
After  ohmic  metal  definition,  the  metal  is  alloyed  to  provide  low  contact 
resistances.  Proton  isolation  implant  follows,  which  is  required  to  increase 
the  resistivity  of  the  substrate  thereby  decreasing  the  effect  of  backgating  in 

LSI  circuits.  The  gate  contact  areas  are  defined  similar  to  the  ohmic  contacts 

and  Ti/Pt  metal  is  deposited  and  lifted  to  form  the  Schottky  contacts. 

D)  A  first  level  SiOxNy  dielectric  is  deposited,  the  wafer  is 
patterned  and  the  metal- 1  pattern  is  recessed  through  the  dielectric  vias  to  the 
contact  metals.  Metal  1  is  then  evaporated  and  lifted.  This  results  in  contact 
pads  for  our  DRPs  and  the  process  was  terminated  at  this  point. 

6.3.3  Automated  Measurements 

DC  parametric  tests  were  performed  on  one  row  of  6400  FETs  on  each 

wafer  of  lots  1  &  2  A  Data  General  MV  4000  computer  was  programmed  to 
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control  the  testing,  analyze  and  store  the  results.  The  Data  General  controlled 
a  HP  41 45 A  Semiconductor  Analyzer  used  in  conjunction  with  an  Eleclroglas 
1034x  automatic  prober,  (see  Fig.  6.3) 


Fig.  6.3.  Automatic  measurement  system  block  diagram. 

Two  scans  of  current  vs  voltage  were  made:  Ijs  vs  Vg  under  constant 
and  Ijs  ''s  under  constant  Vg,  for  each  FET.  This  data  is  digitized  and 
stored  on  the  MV4000  for  later  analysis.  An  analysis  program  performs  a  least 
squares  fit  to  determine  the  parameters  Vth,  Idss,  f^on>  K,  Gm .  Ileak-  Next  a 
plotting  program  graphs  the  resulting  data  vs  their  relative  position  across 
the  3  inch  wafer,  numbering  the  FETS  from  1-6400.  The  data  is  then  averaged 
over  each  DRP  eliminating  any  unwanted  correlation  associated  with  drain 
position  or  number. 

It  will  be  seen  from  the  selenium  implant  data  that  no  significant 
differences  arc  observed  for  Vx  between  the  various  starting  materials. 
However,  the  silicon  implanted  FETs  yielded  data  with  good  correlation  to 
material  properties. 
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6.4  Experimental  Results 

We  have  measured  the  fluctuation  in  MESFET  device  parameters  using  a 
Dense  Row  Pattern  (DRP)  of  MESFETs.  The  devices  have  a  3  pm  gate  length  a 
5nim  gate  width  and  arc  spaced  on  10pm  centers.  This  allows  up  to  6400  FET's  to 
be  measured  on  a  line  across  a  three  inch  wafer.  A  complete  set  of  Process 
Control  Monitors  (PCM)  is  available  for  each  set  of  800  MESFETs.  The  test 
structures  are  fabricated  using  a  simplified,  six  mask  process  with  all  process 
variables  except  material  and  channel  implants  held  constant  to  minimize 
process  induced  defects.  The  MESFETs  arc  fabricated  in  a  foundry  with  a 
proven  process  (Rockwell  MRDC).  This  approach  minimizes  unknown 
variables,  thereby  allowing  enhancement  of  the  effect  of  material  and 
channel  implant  variations  to  be  studied,  while  allowing  statistically 
significant  quantities  of  data  to  be  taken.  Each  wafer  has  6400  MESFETs 
measured  and  characterized  using  an  automatic  tester  and  methods  developed 
previously  at  Rockwell  and  described  above. 

We  have  chosen  to  study  a  matrix  of  material  variation  using  both  Si  and 
Sc  channel  implants.  This  is  to  examine,  in  detail,  their  effect  on  threshold 
voltage  fluctuations  for  the  two  different  channel  implants,  allowing 
elucidation  of  the  effect  described  above.  Material  variations  include  both  low, 
GaAsdn),  dislocation  and  high  dislocation  boulcs  with  and  without  whole  boule 

annealing. 

6,4.1  Correlation  Experiments 

In  Table  6.1  we  show  a  summary  of  the  measurements  completed  to  date. 
Each  value  for  Vx  and  AVj  represents  the  global  average  of  6400  MESFETs 
measured  on  each  wafer.  The  measurement  pattern  is  shown  in  Fig.  6.4.  Each 
reticule  in  Pig.  6.4  consists  of  800  MESFETs  an  a  PCM.  Two  different  substrate 
rnaicrjals  arc  used,  MA7041  a  whole  boule  annealed,  high  dislocation  density 

source  (FT’D  =  .6x10^^  cm'^)  and  R169,  a  Rockwell  grown,  non-whole  boule 
annralcd  low  dislocation  density  source  (EPD  =  10^  cm‘2).  Data  arc  shown  for 
hoih  .Se  and  Si  implanted  channels.  It  is  .seen  from  these  data  that  the  average 
fliK  iu.ilinn  in  A  V  j  is  about  twice  as  large  for  the  Si  implanted  wafers  when 

compared  to  the  Sc  im[)lantcd  ones. 

It  can  al  ()  be  noted  from  Table  6,1  rr,al  both  the  high  and  low  dislocation 
density  matt  ials  exhibit  similar  threshold  fluctuations  when  using  the  same 
ion  for  the  channel  implant.  This  can  be  understood  from  Fig.  6.4  where  we 

have  plotted  the  results  of  Monte-Carlo  calculations  for  AVj  versus  Etch  Pit 
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Density  (EPD)  with  the  base  activation,  Tlo.  as  a  variable.  The  calculations  arc 
for  145  keV,  2.2x10^^  enf^  Si  and  320  keV,  2.2x10^2  cm'^  Se  channel  implants 
through  a  75  nm  Si3N4  encapsulant,  the  conditions  used  in  the  device 
fabrication.  The  interesting  result  of  this  calculation  is  that  for  certain 
conditions,  namely  low  activation  and  high  sensitivity  to  the  presence  of  a 
dislocation  (large  Aq)  a  drifted  gaussian  distribution  results  which  is 

Table  6.1.  Global  average  of  Vx  and  AVj  for  Se  and  Si 

implanted  MESFET's.  Both  whole  boulc  annealed 
high  dislocation  (5x10“^  cm'2,MA  7041)  and  non- 
whole  boulc  annealed  low  dislocation  (10^ 
cm'2,R169)  results  are  listed. 


Boulc 

No. 

Wafer  No. 

Channel 

Implant 

Ion 

Vx 

Se 

DVx 

% 

Vx 

Si 

DVx 

% 

R169 

28 

-0.557 

0.035 

6.28 

In  Doped 

29 

-0.537 

0.033 

6.15 

31 

-0.519 

0.031 

5.97 

32 

-0.534 

0.044 

8.24 

30 

-0.339 

0.066 

16.56 

33 

-0.423 

0.060 

14.03 

.VIA  7041 

148 

-0.556 

0.041 

7.37 

High 

Disk 

150 

-0.455 

0.024 

5.27 

151 

-0.551 

0.028 

5.01 

64 

-0.319 

0.050 

15.64 

65 

-0.424 

0.055 

12.92 

strongly  influenced  by  the  base  activation.  For  the  Si  channel  implants  we 

calculate  a  value  for  r]  0  of  =  64%,  and  for  Sc  =  80%.  Using  average  EPDs  of 
5x10“^  and  10^  cm'^’  it  can  be  seen  from  Fig.  6.5  that  almost  identical  AVjs,  in 
(he  range  of  about  50  mV  result,  in  reasonable  agreement  with  the 

experimental  values.  The  results  of  a  similar  calculation  for  the  Sc  implants 
arc  also  shown  on  this  figure,  however  since  the  base  activation  is  such  that 
little  threshold  fluctuations  arc  seen,  it  is  difficult  to  differentiate  between  the 
fluctuations  caused  by  processing  variables  versus  those  caused  by 
dislocations. 
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Fig.  6.4.  Diagram  of  the  measurement  pattern  and  its  location  with 

respect  to  the  wafer  orientation.  A  total  of  6400  FETs  are 
measured  on  each  wafer.  Each  reticule  contains  800  FETs 
and  one  PCM. 

In  Fig  6.5  we  show  selected  Vj  data  taken  for  both  Se  and  Si  channel 
implants  into  both  high  and  low  EPD  material.  The  higher  fluctuation  in  the  Si 
threshold  voltages  (Fig.  6.5  b  and  c)  are  quite  obvious  in  this  data.  In  Fig.  6.5.d, 

a  negative  step  in  threshold  voltage  is  seen  beginning  at  about  FET  2000  and 

disappears  again  at  about  FET  5000.  This  corresponds  well  with  the  value  of  EL2 
mapped  for  wafers  grown  from  the  same  boule.  In  Fig.  6.6,  we  show  a  plot  of 
EL2  versus  wafer  location  for  an  unannealed  In  doped  Rockwell  wafer  taken 
by  optical  absorption  at  1100  nm.  By  comparing  the  measurement  pattern  of 
Fig.  6.6  with  the  measured  EL2  map,  a  distinct  correlation  can  be  made  between 
the  measured  Vt  and  the  EL2  density,  particularly  at  the  center  of  the  wafer 

which  is  seen  to  be  rich  in  EL2.  Since  EL2  is  generally  believed  to  be  related  to 

the  As  antisite,  AsGa>  ihe  central  region  in  the  GaAs(In)  wafer  should  be  rich 
in  Gava-  T^his  will  increase  the  silicon  activation,  resulting  in  a  more  negative 
Vj  as  seen  in  the  experimental  results. 

6.4.2  Backgating  Measurements 

A  consiuerable  effort  in  programming  of  the  testing  and  analysis  was 
made  on  the  problem  of  backgating.  Whereas  previous  backgating  work 
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Fig.  6.5.  Selected  Vt  data  taken  for  both  Sc  and  Si  channel  implants 

into  both  high  (5x10^^  cm'^)  and  low  (10^  cm"^)  dislocation 
material.  Approximately  6400  devices  are  measured  and 
their  important  parameters  extracted  numerically  for  this 
test. 

reported  in  the  literature  is  done  on  single  devices,  an  attempt  was  made  to 
measure  the  effect  over  6400  FETs  as  was  done  for  the  other  DC  parameters. 
Although  many  test  runs  were  made  and  data  taken,  the  results  arc  not 
promising.  There  arc  two  significant  reasons  for  this:  1)  A  good  model  of  the 
backgating  phenomenon  has  not  been  defined,  that  is.  to  determine  Vj,  for 
example  a  least  squares  fit  of  the  square  root  of  Ids  ''s  Vgs  is  made  on  digitized 
data.  No  similar  standard  exists  for  backgating.  Some  investigators  report  a 
10%  shift  in  Vj ,  and  some  a  10%  shift  in  Idss-  Likewise  there  is  no  standard 
geometry  for  the  placement  of  the  backgalc.  2)  The  initial  mask  set,  while 
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Fig.  6.6.  EL2  map  of  an  In  doped,  low  dislocation  wafer  from  the 

same  boule  as  the  device  wafers.  A  region  of  high  EL2  (and 
dislocation  density)  is  seen  in  the  center  of  this  wafer. 
Increasing  darkness  in  the  gray  scale  corresponds  to 
increasing  EL2  content. 

concentrating  on  minimizing  processing  risks,  was  made  using  only  one 
metalization  layer.  This  forces  the  interconnect  metal  to  lye  on  bare  GaAs.  In 
the  backgating  measurements  this  affected  the  data  by  causing  the  results  to 
be  dependent  on  the  drain  number  of  the  FET  in  the  DRP.  There  remains  a 
need  in  the  GaAs  industry  for  a  more  comprehensive  backgate  model  and 
testing  procedures. 

6.4.3  Etching  Experiments 

Four  of  the  wafers  in  lot  #2  were  to  be  used  to  correlate  the  dislocation 
density  with  the  measured  electrical  results.  First  the  interconnect  metals 
were  removed  by  ion  milling.  This  left  a  demarcation  on  the  bare  GaAs 

substrate  of  the  gate  location  of  each  FET.  These  wafers  were  then  immersed  in 
fresh  molten  KOH  at  450C  for  20-25  minutes.  This  is  a  standard  etchant 
commonly  used  to  reveal  dislocations.  The  wafers  were  then  cooled  and  rinsed 
in  deionized  water.  Photomicrographs  were  taken  of  an  entire  row  containing 
the  electrically  tested  FETs.  Three  sets  of  photos  were  analyzed,  R169-32&33 
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(In-low  dislocation),  and  MA7041-64  (standard  LEC).  A  DRP  consists  of  units  of 
10  FETs  with  common  source  and  gate.  The  parameters  of  each  set  of  10  FETs 
were  averaged  along  the  DRP.  This  set  of  640  DRPs  are  compared  to  an 
estimate  made  for  the  dislocation  density.  The  estimate  of  the  dislocation 
density  was  taken  as  follows;  the  number  cf  dislocations  within  a  150  m  circle 
centered  about  each  DRP  subset  was  counted  for  the  640  DRPs.  Due  to  the  high 
number  of  dislocations  in  standard  material,  only  three  fields  were  counted  on 
MA7041-64.  This  estimate  of  the  dislocation  density  was  plotted  vs  the  standard 
deviation  of  the  threshold  voltage  of  each  of  the  FET  subsets.  Figure  6.7  shows 
the  results  of  these  plots.  Figure  6.7  shows  the  standard  deviation  of 
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Fig.  6.7.  Measured  standard  deviation  of  VT  and  EPD  for  wafer  R169-32. 


Vx  across  wafer  R1 69-32  and  shows  the  dislocation  density  across  these  same 
FETs.  Note  where  there  is  a  jump  in  dislocation  density  there  is  a  substantial 
dip  (more  negative)  in  the  corresponding  Vj.  Figure  6.8  shows  a  photo  of  that 
region.  This  represents  the  only  region  in  the  low  dislocation  material 
where  the  dislocation  cells  arc  near  the  the  measured  FETs.  This  cellular 
structure  observed  is  due  to  constitutional  supercooling  of  the  In-GaAs  alloy. 
That  is,  due  to  the  segregation  coefficient  of  In  in  GaAs,  as  the  bulk  growth 
proceeds  In  is  rejected  from  the  growth  interface  resulting  in  an  increase  in 
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Fig.  6.8.  Photographs  of  etched  Wafer  R169-32.  The  balck  dots  represent 
etch  pits.  The  outline  of  the  devices  can  be  seen. 


Fig,  6.9.  Reflection  X-ray  topograph  of  In  doped  Wafer  R169-18.  The 
disrupted  area  in  the  low-r  right  indicates  a  loss  of  crystallinity. 
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crystalline  structure  results.  This  phenomenon  can  be  observed  in  Figure  6.9, 
a  reflection  X-ray  photograph  of  wafer  R169-18.  The  disrupted  area  on  the  X- 
ray  photo  is  an  indication  of  the  loss  of  crystallinity  in  that  area.  This  wafer  is 
from  the  same  boule  as  R169-32  only  it  is  14  wafers  toward  the  seed  end  in 
front  of  it.  This  breakdown  induces  dislocations  which  trail  into  the  line  of 
FETs  as  observed  in  Fig  6.8  above.  The  FETs  measured  electrically  were  not  in 
the  area  of  the  breakdown. 

6.5  Conclusions 

The  results  we  have  ob.served  indicate  that  a  strong  correlation  of  Vj 
with  dislocation  density  does  exist.  However,  the  number  of  process  and 
material  variables  which  are  uncontrolled  tend  to  mask  the  correlation.  A  way 
to  resolve  these  cross  correlated  results  is  to  increase  the  amount  of  data  for 
statistical  purposes.  The  model  developed  during  this  effort  has  predicted  the 
results  quite  adequately.  Processing  of  Lot  #3  would  certainly  strengthen  the 
final  conclusions,  as  it  would  include  both  WHOBA  In  doped  material  (to 
redistribute  the  EL2)  and  VGF  AT&T  material  (low  EPD  uniform  EL2). 
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